Classifying quantum phases using MPS and PEPS 
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We give a classification of gapped quantum phases of one-dimensional systems in the framework of 
Matrix Product States (MPS) and their associated parent Hamiltonians, for systems with unique as 
well as degenerate ground states, and both in the absence and presence of symmetries. We find that 
without symmetries, all systems are in the same phase, up to accidental ground state degeneracies. 
If symmetries are imposed, phases without symmetry breaking (i.e., with unique ground states) 
are classified by the cohomology classes of the symmetry group, this is, the equivalence classes of 
its projective representations, a result first derived in [X. Chen, Z.-C. Gu, and X.-G. Wen, Phys. 
Rev. B 83, 035107 (2011): larXiv: 1008.3745] . For phases with symmetry breaking (i.e., degenerate 
ground states), we find that the symmetry consists of two parts, one of which acts by permuting 
the ground states, while the other acts on individual ground states, and phases are labelled by 
both the permutation action of the former and the cohomology class of the latter. Using Projected 
Entangled Pair States (PEPS), we subsequently extend our framework to the classification of two- 
dimensional phases in the neighborhood of a number of important cases, in particular systems with 
unique ground states, degenerate ground states with a local order parameter, and topological order. 
We also show that in two dimensions, imposing symmetries does not constrain the phase diagram 
in the same way it does in one dimension. As a central tool, we introduce the isometric form, a 
normal form for MPS and PEPS which is a renormalization fixed point. Transforming a state to its 
isometric form does not change the phase, and thus, we can focus on to the classification of isometric 
forms. 



I. INTRODUCTION 
A. Background 

Understanding the phase diagram of correlated quan- 
tum many-body systems, that is, the different types of 
order such systems can exhibit, is one of the most impor- 
tant and challenging tasks on the way to a comprehen- 
sive theory of quantum many-body systems. Compared 
to classical statistical models, quantum systems exhibit a 
much more complex behavior, such as phases with topo- 
logical order as in the fractional quantum Hall effect, 
which cannot be described using Landau's paradigm of 
local symmetry breaking. 

In the last years, tensor network based ansatzes, such 
as Matrix Product States^ (MPS) and Projected Entan- 
gled Pair States^ (PEPS), have proven increasingly suc- 
cessful in describing ground states of quantum many- 
body systems. In particular, it has been shown that 
MPS and PEPS can approximate ground states of gapped 
quantum systems efficiently^— this is, not only the 
ground state of systems with local order, but also, for 
instance, the ground states of topological insulators are 
well represented by those states. Since MPS and PEPS 
provide a characterization of quantum many-body states 
from a local description, they are promising candidates 
for a generalization of Landau's theory. Moreover, they 
can be used to construct exactly solvable models, as ev- 
ery MPS and PEPS appears as the exact ground state of 
an associated parent Hamiltonianii^r— 

In this paper, we apply the framework of MPS and 
PEPS to the classification of gapped quantum phases 



by studying systems with exact MPS and PEPS ground 
states. Here, we define two gapped systems to be in the 
same phase if and only if they can be connected by a 
smooth path of gapped local Hamiltonians. Along such 
a path, all physical properties of the state will change 
smoothly, and as the system follows a quasi-local evolu- 
tion^ global properties are preserved. A vanishing gap, 
on the other hand, will usually imply a discontinuous be- 
havior of the ground state and affect global properties of 
the system. In addition, one can impose symmetries on 
the Hamiltonian along the path, which in turn leads to a 
more refined classification of phases. 

In the presence of symmetries, the above definition of 
gapped quantum phases can be naturally generalized to 
systems with symmetry breaking, i.e., degenerate ground 
states, as long as they exhibit a gap above the ground 
state subspace. Again, ground states of such systems 
are well approximated by MPS and PEPS, which justi- 
fies why we study those phases by considering systems 
whose ground state subspace is spanned by MPS and 
PEPS. On the other hand, the approach will not work 
for gapless phases, as the MPS description typically can- 
not be applied to them. 



B. Results 

Using our framework, we obtain the following classifi- 
cation of quantum phases: 

For one-dimensional (ID) gapped systems with a 
unique ground state, we find that there is only a sin- 
gle phase, represented by the product state. Imposing a 
constrain in form of a local symmetry U g (with symme- 
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try group G 9 g) on the Hamiltonian leads to a more 
rich phase diagram. It can be understood from the way 
in which the symmetry acts on the virtual level of the 
MPS, U g = V g ® V g , where the Vg are projective repre- 
sentations of G. In particular, different phases under a 
symmetry are labelled by the different equivalence classes 
of projective representations V g of the symmetry group 
GfiS which are in one-to-one correspondence to the ele- 
ments of its second cohomology group H 2 (G,U(1)) (this 
has been previously studied in Ref. fill). 

For one-dimensional gapped systems with degenerate 
ground states, we find that in the absence of symme- 
tries, all systems with the same ground state degeneracy 
can be transformed into another along a gapped adia- 
batic path. In order to make these degeneracies stable 
against perturbations, symmetries need to be imposed 
on the Hamiltonian. We find that any such symmetry 
decomposes into two parts, P g and Wh- Here, P g acts 
by permuting the symmetry broken ground states of the 
system. To describe Wh, choose a "reference" ground 
state, and let H C G be the subgroup for which Ph acts 
trivially on the reference state: Then, Wh (h G H) is 
a unitary symmetry of the reference state, which again 
acts on the virtual level as Wh — Vh ® Vh, with Vh a pro- 
jective representation of H . Together, P g and Wh form 
an induced representation. The different phases of the 
system are then labelled both by the permutation action 
P g of U g on the symmetry broken ground states (or alter- 
natively by the subgroup H C G for which Ph leaves the 
reference state invariant), and by the equivalence classes 
of projective representations Vh of H. 

Our classification of phases is robust with respect to 
the definition of the gap: Two systems which are within 
the same phase can be connected by a path of Hamilto- 
nians which is gapped even in the thermodynamic limit; 
conversely, along any path interpolating between systems 
in different phases the gap closes already for any (large 
enough) finite chain. On the other hand, we demonstrate 
that the classification of phases is very sensitive to the 
way in which the symmetry constraints are imposed, and 
we present various alternative definitions, some of which 
yield more fine-grained classifications, while others result 
in the same classification as without symmetries. In par- 
ticular, we also find that phases under symmetries are 
not stable under taking multiple copies, and thus should 
not be regarded a resource in the quantum information 
sense. 

Parts of our results can be generalized to two dimen- 
sions (2D), with the limitation that we can only prove 
gaps of parent Hamiltonians associated with PEPS in 
restricted regions. These regions include systems with 
unique ground states and with local symmetry breaking 
as well as topological models. We show that within those 
regions, these models label different quantum phases, 
with the product state, Ising Hamiltonians, and Kitaev's 
double models^ as their representatives. We also find 
that in these regions, imposing symmetries on the Hamil- 
tonian does not alter the phase diagram, and, more gen- 



erally, that symmetry constraints on two- and higher- 
dimensional systems must arise from a different mecha- 
nism than in one dimension. 

As a main tool for our proofs, we introduce a new stan- 
dard form for MPS and PEPS which we call the isometric 
form. Isometric forms are renormalization fixed points 
which capture the relevant features of the quantum state 
under consideration, both for MPS and for the relevant 
classes of PEPS. Parent Hamiltonians of MPS can be 
transformed into their isometric form along a gapped 
path, which provides a way to renormalize the system 
without actually blocking sites. This reduces the classi- 
fication of quantum phases to the classification of phases 
for isometric MPS/PEPS and their parent Hamiltonians, 
which is considerably easier to carry out due to its addi- 
tional structure. 



C. Structure of the paper 

The paper is structured as follows: In Section II, we 
prove the results for the one-dimensional case: We start 
by introducing MPS (Sec. II A) and parent Hamilto- 
nians (Sec. II B), and define phases without and with 
symmetries (Sec. II C). We then introduce the isomet- 
ric form and show that the problem of classifying phases 
can be reduced to classifying isometric forms (Sec. II D). 
Subsequently, we first classify ID phases without sym- 
metries (Sec. II E), then phases of systems with unique 
ground states under symmetries (Sec. II F), and fi- 
nally phases of symmetry broken systems under symme- 
tries (Sec. II G). 

In Section III, we discuss the 2D scenario. We start 
by introducing PEPS (Sec. Ill A) and characterize the 
region in which we can prove a gap (Sec. Ill B). We then 
classify PEPS without symmetries (Sec. Ill C) and show 
that symmetries don't have an effect comparable to one 
dimension (Sec. Ill D). 

Section IV contains discussions of various topics which 
have been omitted from the preceding sections. Most im- 
portantly, in Sec. IV B and C, we discuss various ways 
in which phases, in particular in the presence of symme- 
tries, can be defined, and the way in which this affects 
the classification of phases, and in Sec. IIVD1 we provide 
examples illustrating our classification. 



II. RESULTS IN ONE DIMENSION 

In this section, we will derive the classification of 
phases for one-dimensional systems both with unique and 
with degenerate ground states, and both in the absence 
and the presence of symmetries. We will start by giv- 
ing the necessary definitions — we will introduce Matrix 
Product States and their parent Hamiltonians, and define 
what we mean by phases both without and with symme- 
tries. Then, we will state and prove the classification 
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FIG. 1. MPS are constructed by applying a linear map V to 
maximally entangled pairs \ojd) := Yst=i l*>*) °f bond dimen- 
sion D. 

of phases for the various scenarios (with some technical 
parts placed in appendices). 

Note that for clarity, we will keep discussions in this 
section to a minimum. Extensive discussion of various 
aspects (motivation of the definitions, alternative defini- 
tions, etc.) can be found in Sec. IIVI 



A. Matrix Product States 

In the following, we will study translational invari- 
ant systems on a finite chain of length TV with periodic 
boundary conditions. While we do not consider the ther- 
modynamic limit, we require relevant properties such as 
spectral gaps to be uniform in N. 



1. Definition of MPS 

Consider a spin chain (<C d )® N . A (translational invari- 
ant) Matrix Product State (MPS) \n\P\) of bond dimen- 
sion D on (<C d )® N is constructed by placing maximally 
entangled pairs 

D 

\lo d ) := M 

i=l 

between adjacent sites and applying a linear map V : 
C D (g> C D -> C d , as depicted in Fig. [1^ this is, \fi[V]) = 
V® N \uo D )® N . The map V is sometimes called the MPS 
projector (though it need not be a projector). 



2. Standard form 

The definition of MPS is robust under blocking sites: 
Blocking k sites into one super-site of dimension d k gives 
a new MPS with the same bond dimension, with projec- 
tor V' = r P® k \uiD)® { - k ~ 1 \ where the maximally entangled 
pairs are again placed as in Fig. [TJ 

By blocking and using gauge degrees of freedom (in- 
cluding changing £>), any MPS which is well-defined in 
the thermodynamic limit can be brought into its stan- 
dard form, where V is supported on a "block-diagonal" 
space, (kcr V) 1 - = Hi 8 • • • © H A , with 

U a = span{|i, j) : Ca-i < i,j < (a} • (1) 



Here, = Co < Ci < " ' ' < Ca = D gives a partitioning 
of 1, ... ,D; we define Dj := Q — Ci-i- The case where 
A = 1 (i.e. V is injective) is called the injective case, 
whereas the case with A. > 1 is called non- injective. 

In the following, we will assume w.l.o.g. that all MPS 
are in their standard form — note that this might involve 
the blocking of sites of the original model. Moreover, we 
will impose that V is surjective, this can be achieved by 
restricting C d to the image of V. 

B. Parent Hamiltonians 

For any MPS (in its standard form) one can construct 
local parent Hamiltonians which have this MPS as their 
ground state. A (translational invariant) parent Hamil- 
tonian 

N 

H = J2h(i,i + l) 

consists of local terms h(i, i + 1) = h > acting on two 
adjacent sites (i, i + 1) whose kernels exactly support the 
two-site reduced density operator of the corresponding 
MPS, i.e., 

kerh = {V®V){C D ® \u D )®C D ) . (2) 

By construction, H > and H\fi[P]) = 0, i.e., l^tfP]) is 
a ground state of H. It can be shown that the ground 
state space of H is „4-fold degenerate, and spanned by 
the states |ai['P|'h ( J)i where V\n a is the restriction of V 
to l-Lar&L This associates a family of parent Hamilto- 
nians to every MPS, and by choosing h a projector, the 
mapping between MPS and Hamiltonians becomes one- 
to-one. 

Given this duality between MPS and their parent 
Hamiltonians, we will use the notion of MPS and par- 
ent Hamiltonians interchangingly whenever appropriate. 

C. Definition of quantum phases 

Vaguely speaking, we will define two systems to be in 
the same phase if they can be connected along a con- 
tinuous path of gapped local Hamiltonians, possibly pre- 
serving certain symmetries; here, gapped means that the 
Hamiltonian keeps its spectral gap even in the thermo- 
dynamic limit. The intuition is that along any gapped 
path, the expectation value of any local observable will 
change smoothly, and the converse is widely assumed to 
also holdJ^ 

The rigorous definitions are as follows. 

1. Phases without symmetries 

Let Hq and Hi be a family of translational invari- 
ant gapped local Hamiltonians on a ring (with periodic 
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boundary conditions). Then, we say that Hq and Hi are 
in the same phase if and only if there exists a finite k such 
that after blocking k sites, Hq and Hi are two-local, 

N 

H v = ^h p {i,i + l) , p = 0,l, 

i=l 

and there exists a translational invariant path 

N 

H i =£mm+i) , o<7<i, 

i=l 

with two-local /i 7 such that 

i) ft-o = ^7=0 j fii = h 1= i 

ii) II^IIop < 1 

iii) hy is a continuous function of 7 

iv) i? 7 has a spectral gap above the ground state mani- 
fold which is bounded below by some constant A > 
which is independent of TV and 7. 

In other words, two Hamiltonians arc in the same phase 
if they can be connected by a local, bounded-strength, 
continuous, and gapped path. 

Note that this definition applies both to Hamiltonians 
with unique and with degenerate ground states. 

2. Phases with symmetries 

Let H p (p = 0, 1) be a Hamiltonian acting on Hf N , 
Hp = C dp , and let UJ} be a linear unitary representation 
of some group G 3 g on H p . We then say that U g is 
a symmetry of H p if [H p , (Ug 3 )®"} = for all g g G; 
note that is only defined up to a one-dimensional 
representation of G, f/| <-> e^aU?. Then, we say that 
Hq and Hi are in the same phase under the symmetry 
G if there exists a phase gauge for Ug and Ug and a 
representation U = Ug © Ug © C/P ath of G on a Hilbcrt 
space "H = 'Hoffi'Hiffi'Hpath, and an interpolating path H 1 
on H with the properties given in the preceding section, 
such that [H 1: Uf N ] = 0, and where Hq and Hi are 
supported on Hq and Hi, respectively. 

There are a few points to note about this definition: 
First, we allow for an arbitrary representation of the sym- 
metry group along the path; we will discuss in Sec. IIV CI 
why this is not a restriction. Second, we impose that 
Hq and Hi are supported on orthogonal Hilbcrt spaces: 
This allows us to compare e.g. the spin-I AKLT state 
with the spin-0 state under SO (3) symmetry, but we will 
impose this even if the two representations are the same; 
we will discuss how to circumvent this in Sec. IIV Cl Note 
that, just as without symmetries, this definition should 
be understood after an appropriate blocking of sites. 

3. Robust definition of phases 

In addition to the properties listed in the definition of 
phases in Sec. Ill C II above, one usually requires a phase 



to be robust, i.e., an open set in the space of allowed 
Hamiltonians: For every Hamiltonian 

N 

H = J2 h (h i + 1 ) 

i=l 

there should be an e > such that 

N 

H = !) + f-k(i,i + 1) 

i=i 

is in the same phase for any bounded-strength k(i,i + 1) 
which obeys the required symmetries. 

We are not going to rigorously address robustness 
of phases in the present paper; however, it should be 
pointed out that, in the absence of symmetries, Hamil- 
tonians with degenerate MPS ground states do not sat- 
isfy this property: In its standard form, the different 
ground states of a non-injective MPS are locally sup- 
ported on linearly independent subspaccs, and we can 
use a translational invariant local perturbation e k(i) to 
change the energy of any of the ground states propor- 
tionally to eN, thereby closing the gap for a N oc I/e. 
On the other hand, in the presence of a symmetry which 
permutes the different ground states (such a symmetry 
always exists), those perturbations are forbidden, and 
the phase becomes stable. (A rigorous stability proof for 
MPS phases would make use of the stability condition 
for frustration-free Hamiltonians proven in Ref. [IH and 
its generalization to symmetry-broken phases analogous 
to the one discussed in Ref. [LQ The condition is trivially 
satisfied by the renormalization fixed points of MPSjii 
and using the exponential convergence of MPS to their 
fixed point ji£ the validity of the stability condition, and 
thus the stability of MPS phases, follows.) 

Therefore, when classifying phases of systems with de- 
generate ground states, one should keep in mind that in 
order to make this a robust definition, a symmetry which 
protects the degeneracy is required. 



4- Restriction to parent Hamiltonians 

We want to classify the quantum phases of gapped 
Hamiltonians which have exact MPS ground states (or, 
in the case of degeneracies, the same ground state sub- 
space as the corresponding parent Hamiltonian). For- 
tunately, with our definition of phases it is sufficient to 
classify the phases for parent Hamiltonians themselves: 
Given any two gapped Hamiltonians H and H' which 
have the same ground state subspace, the interpolating 
path jH + (1 — j)H' has all desired properties, and in 
particular it is gapped. Note that this also shows that 
all parent Hamiltonians for a given MPS are interchange- 
able. 
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FIG. 2. Construction of the interpolating path for MPS and 
parent Hamiltonians. Instead of interpolating between the 
MPS ) an d \ipi) directly (dotted line), we first show how 
to interpolate each of the two states towards a standard from 
\ipp), the isometric form, and then construct an interpolating 
path between the isometric forms. Note that using the parent 
Hamiltonian formalism, any such path in the space of MPS 
yields a path in the space of Hamiltonians right away. 



D. The isometric form 

1. Reduction to a standard form 

Given two MPS \fJ,[P p ]), p = 0, 1, together with their 
nearest- neighbor parent Hamiltonians H p , we want to see 
whether H$ and Hi arc in the same phase, i.e., whether 
we can construct a gapped interpolating path. We will do 
so by interpolating between Vq and V\ along a path V 1 , 
in such a way that it yields a path H 1 in the space of par- 
ent Hamiltonians which satisfies the necessary continuity 
and gappedness requirements. 

In order to facilitate this task, we will proceed in two 
steps, as illustrated in Fig. H In a first step, we will 
introduce a standard form for each MPS — the isometric 
form — which is always in the same phase as the MPS 
itself. This will reduce the task of classifying phases to 
the classification of phases for isometric MPS, which we 
will pursue subsequently. 



2. The isometric form 

Let us now introduce the isometric form of an MPS. 
The isometric form captures the essential entanglement 
and long-range properties of the state and forms a fixed 
point of a renormalization procedure,— and every MPS 
can be brought into its isometric form by stochastic local 
operations.— Most importantly, as we will show there 
exists a gapped path in the space of parent Hamiltonians 
which interpolates between any MPS and its isometric 
form. 

Given an MPS |/x["P]), decompose 



V 



Qt 



w 



FIG. 3. Isometric form of an MPS. a) The MPS projector V 
can be decomposed into a positive map Q and an isometric 
map W . b) By removing Q, one obtains the isometric form W 
of the MPS. c) Interpolation to the isometric form is possible 
by letting Q 7 = 7 Q + (1 - j)t. 



define an interpolating path l/ifRJ), V 1 = QjW, with 

Q 7 = 7 Q +(1- 7 )1 , 1> 7 >0; 

note that the path can be seen as a stochastic deforma- 
tion \(J,[P-y}) = Q® N \fi[V ]), cf. Fig. [3] Throughout the 
path, the MPS stays in standard form, and in the non- 
injective case, the blocking pattern which is encoded in 
the structure of ker V stays unchanged. 

Let us now see that l/xp 7 ]) = |/x['Pi]) and its isometric 
form |/x[W]) = | ['Po] ) ar e m the same phase. To this 
end, consider the parent Hamiltonian Hq = ^2 ho(i, i + l) 
of the isometric MPS | [T^o] ) : with Hq a projector. Let 



A, 



and define the "/-deformed Hamiltonian 



H~ 



y h 7 (i, i + 1) , with /i 7 



A 7 /i A 7 > 



Since ho\fi[Po\) — 0, it follows that /i 7 |/x[7 7 7 ]) = (and in 
fact, the kernel of h is always equal to the support of the 
two-site reduced state), i.e., i? 7 is a parent Hamiltonian 
of | M [P-y ] ) . Note that the family H 1 of Hamiltonians is 
continuous in 7 by construction. 

It remains to show that the path i? 7 is uniformly 
gapped, i.e., there is a A > which lower bounds the 
gap of -ff 7 uniformly in 7 and the systems size N: This 
establishes that the l/xp^]) are all in the same phase. 
The derivation is based on a result of Nachtergaele^ (ex- 
tending the results of Ref . [l| for the injective case) , where 
a lower bound (uniform in N) on the gap of parent Hamil- 
tonians is derived, and can be found in Appendix A. The 
central point is that the bound on the gap depends on 
the correlation length £ and the gap of H 1 restricted to 
£ sites, and since both depend smoothly on 7, and £ — > 
as 7 — > 0, a uniform lower bound on the gap follows; for 
the non-injective case, one additionally needs that the 
overlap of different ground states goes to zero as 7 — > 0. 



V = Q W , (3) 

with W an isometry WW^ = 11 and Q > 0, by virtue 
of a polar decomposition of V\ (k er w.l.o.g., we can 
assume < Q < 1 by rescaling V . The isometric form of 
of \fi[V]) is now defined to be \fi[W]), the MPS described 
by W, the isometric part of the tensor V. 

To see that |/x["P]) and are in the same phase, 



3. Isometric form and symmetries 

An important point about the isometric form |/i[Po]) 
of an MPS |/Li[Pi]) is that both remain in the same phase 
even if symmetries are imposed. The reason is that using 
gauge transformations 

V <— »■ V(Y ®{Y~ 1 ) T ) 
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h D 

*-C 

i%H ^#91 ##; 

ka> I^b> 

FIG. 4. The isometric form of an injective MPS consists of 
maximally entangled states \uid) with bond dimension D be- 
tween adjacent sites, where each site is a D 2 -level system. 
The local Hamiltonian terms hP , Eq. Q , acts only on the de- 
grees of freedom corresponding to the central entangled pair. 

— such transformations do not change |/Lt[P]) — any V can 
be brought into a standard form where tr loft ['P' l '7 : '] = 
1 right [cf. Refs. [l] and 0j; in this standard form, any 
symmetry Uf N of the MPS l/^T 3 ]) can be understood as 
some unitary X g acting on the virtual level, 7 

u g v = rx g . 

In the polar decomposition V = QW, Eq. ((3j> , we have 
that Q = y/W\ and thus 

q2 = vv \ = v t p X g X\V^ Ug 

= U\VV^XJ g = UlQ 2 U g . 

This is, for any g £ G the matrices Q 2 and U g are di- 
agonal in a joint basis, and therefore [Q, U g ] = 0, and 
it follows that both the interpolating path l/xp^]) = 
Q® N \(i[Po]) and its parent Hamiltonian H 1 are invari- 
ant under U g . 

E. Phase diagram without symmetries 

The preceding discussion show thats in order to classify 
quantum phases (both without and with symmetries), 
it is sufficient to consider isometric MPS and their par- 
ent Hamiltonians. In the following, we will carry out 
this classification for the scenario where no symmetries 
are imposed. We will find that without symmetries, the 
phase of the system only depends on the ground state de- 
generacy, as any two systems with the same ground state 
degeneracy are in the same phase (and clearly the ground 
state degeneracy cannot be changed without closing the 
gap). 

1. Unique ground state 

Let us start with the injective case where the Hamilto- 
nian has a unique ground state. In that case, the isom- 
etry W is a unitary. We can now continuously undo 
the rotation W — this clearly is a smooth gapped path 
and does not change the phase. This yields the state 
\fi[P = 1]) = \wd)® n , consisting of maximally entan- 
gled pairs of dimension D between adjacent sites which 
can only differ in their bond dimension D, cf. Fig. 2] The 



parent Hamiltonian is a sum of commuting projectors of 
the form 

h D = t - \uj d )(lj d \ (4) 

and thus gapped. (Here, hP should be understood as 
acting trivially on the leftmost and rightmost ancillary 
particle, and non-trivially on the middle two, cf. Fig. [U) 
Moreover, for different D and D' one can interpolate be- 
tween these states via a path of commuting Hamiltonian 
with local terms 

h e = l-\ U (0)MO)\ , (5) 

where \u>(9)) = V9\ujd) + \/l - 8\u) D >). It follows that 
any two isometric injective MPS, and thus any two injec- 
tive MPS, are in the same phase. In particular, this phase 
contains the product state as its canonical representant. 

Note that in the previous derivation, the dimensions 
of the local Hilbert spaces differ. One can resolve this by 
thinking of both systems as being embedded in a larger 
Hilbert space, or by adding ancillas; we discuss this fur- 
ther is Sec. [iVl 

2. Degenerate ground states 

Let us now consider the case of of non-injective MPS, 
i.e., systems with degenerate ground states? First, con- 
sider the case with block sizes D a := dimHa = 1, i.e., 
V = Vgwl = J2 a l a X a i a \ ( U P to a rotation of the phys- 
ical system): This describes an .4- fold degenerate GHZ 
state with commuting Hamiltonian terms 

h = 1 — \a, a)(a, a\ . (6) 

a 

For D a ^ 1, we have that (again up to local rotations) 
V = l a X a ' a l ® t Dl J 

Oi 

where acts on T-L a in ([1]). This describes a state 

^| a ,...,a)®|^J 8JV (7) 

a 

(cf. Fig. [5]), i.e., a GHZ state with additional local en- 
tanglement between adjacent sites, where the amount of 
local entanglement D a can depend on the value a of the 
GHZ state. [Since the D a can be different, the correct 
interpretation of is \u Dl )® N © • • • © \uj Da )® n , with 
the | a) labelling the direct sum components locally.] The 
corresponding Hamiltonian is a sum of the GHZ Hamil- 
tonian (|6|) and terms 

^2\a)(a\»(l-\tJ Da )(u Da \) , 

a 

which arc responsible for the local entanglement. This 
Hamiltonian commutes with the GHZ part (this can be 
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V \a,...,a) 

|wdJ\,„* |WD Q )*,.' 

FIG. 5. Structure of the isometric form for non-injective 
MPS: The state consists of a GHZ state ^ Q |a, . . . , a) (gray) 
and maximally entangled pairs between adjacent sites \oJr> a ), 
where the bond dimension D a can couple to the value a of 
the GHZ state. 

understood by the fact that the GHZ state breaks the 
local symmetry, i.e., there is a preferred local basis), and 
one can again interpolate to the pure GHZ state with 
D' a = 1 analogously to the injective case, Eq. ©. 



3. Summary: Classification of phases 

Together, we obtain the following result on the clas- 
sification of ID phases in the absence of symmetries: 
Any two systems with MPS ground states with the same 
ground state degeneracy are in the same phase; the 
canonical representant of these phases are the product 
and the GHZ-type states, respectively. Since the ground 
state degeneracy cannot be smoothly changed without 
closing the gap, this completes the classification. 



F. Phase under symmetries: unique ground state 

Let us now discuss the classification of phases in the 
presence of symmetries, as denned in Sec. Ill C 21 We will 
in the following first discuss the case of injective MPS, 
i.e., systems with unique ground states. Note that we 
will consider the symmetry of the blocked MPS, but as 
we have argued when defining phases under symmetries, 
this does not affect the classification. 

An important prerequisite for the subsequent discus- 
sion is the observation that any MPS has a gauge degree 
of freedom 

e i °\n[V]) = \ t i[e i *V{Y®Y*)}) , (8) 

where Y is right-invertible (i.e., there exists Y" 1 s.th. 
YY^ 1 = Id, but Y need not be square), and Y* = 
(y _1 ) T . Conversely, it turns out that any two V repre- 
senting the same state are related by a gauge transfor- 
mation ||5J|. 



1. Projective representations and the classification of 
phases 

Let U g be a linear unitary representation of a sym- 
metry group G. We start from the fact that for any 
t/ 9 -invariant MPS l/^T^]) and parent Hamiltonian, there 



is a standard form for V and a phase gauge for U g such 
that 

v g v = v(y B ®v g ) , (9) 

where the bar denotes the complex conjugate. 7 Here, the 
V g form a projective unitary representation of group G, 
i.e., V g Vh = e iuj ^V gh ; cf. Appendix H for details. As 
V g appears together with its complex conjugate in @, it 
is only defined up to a phase, V g O e ZXg V g , and thus, 
uj(g,h) is only determined up to the equivalence relation 

u(g, h) ~ ui(g, h) + x g h ~ Xg ~ Xh mod 2ir . 

The equivalence classes induced by this relation form a 
group under addition (this is, tensor products of repre- 
sentations), which is isomorphic to the second cohomol- 
ogy group H 2 (G, U(l)) of G over U(l); thus, we will also 
call them cohomology classes. 

In the following, we will show that in the presence of 
symmetries, the different phases are exactly labelled by 
the cohomology class of the virtual realization V g of the 
symmetry U g determined by Eq. Q; this result was pre- 
viously found in Ref. [Til . 

2. Equality of phases 

Let us first show that MPS with the same cohomology 
class for V g in © can be connected by a gapped path. 
We will do so by considering the isometric point, where 
V is unitary; recall that the transformation to the iso- 
metric point commutes with the symmetry, so that @ 
still holds. Then, © can be rephrased as 

Ug = V^UgV = Vg®Vg, 

this is, the action of the symmetry can be understood as 
U g = V g ® V g acting on the virtual system, in a basis 
characterized by V . 

Now consider two MPS with isometric forms Vq and 
Vi, and symmetries 

U^ = VtU° g V = Vg" 8 Vg , 

U] = V\V\Vi = v x g ® V g , 

where V g and V g are in the same cohomology class. We 
can now interpolate between the two MPS with bond 
dimensions Dq and D\ (in the "convenient basis" corre- 
sponding to Ug and U g ) along the path |w(7)) l8Ar , where 

M7)H(l-7)£M>+7 E M>< ( 10 ) 

i=l i=D +l 

which is an MPS with bond dimension Do + D±. Again, 
the parent Hamiltonian along this path is commuting and 
thus gapped, and changes smoothly with 7. This path 
can be understood as a path with symmetry 

(Vg © Vg^ 8 (Vg Q = ^ Q ^ ^ , 
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(cf. Ref.[Iit), with 

£path = (yO g, yl) e (yl g, ^0) 

(this is where equality of the cohomology classes is re- 
quired, since only then {/P ath forms a linear representa- 
tion) . 

3. Separation of phases 

As we have seen, MPS for which V g in Eq. ^ is in the 
same cohomology class fall into the same phase. Let us 
now show that conversely, states with different cohomol- 
ogy classes fall into different phases. We will prove this 
again in the framework of MPS, i.e., we will show that 
there cannot be a smooth MPS path connecting two such 
states. Note that it is clear that the interpolation given 
above cannot work, as now V® © Vg does not form a 
representation any more. 

The idea of the impossibility proof is to consider a 
chain of arbitrary length N and show that along any 
well-behaved path H 7 , V 1 needs to change continuously, 
which results in a continuous change in the way the sym- 
metry acts on the virtual system. In turn, such a contin- 
uous change cannot change the cohomology class. While 
this argument is based on the fact that the chain is finite 
(as the continuity bounds depend on N), it works for ar- 
bitrary system size N; also, our argument implies that in 
order to interpolate between two systems with different 
cohomology classes, i.e., in different phases, the gap of 
the Hamiltonian will have to close for a finite chain, and 
not only in the thermodynamic limit. (This can be un- 
derstood from the fact that along an MPS path, the vir- 
tual representation of the symmetry is well defined even 
for finite chains, so that it cannot change without closing 
the gap. While we believe that cohomology classes la- 
bel gapped phases beyond MPS, this will likely not hold 
exactly for finite chains, thus leaving the possibility of a 
higher order phase transition when interpolating beyond 
MPS.) We will now proceed by fixing some 7 along the 
path and show the continuity in an environment 7 + cky. 

Fix some 7, with corresponding ground state |m["P 7 ]) 
of H-f, where P 7 is w.l.o.g. in its standard form where 
^JgPi = ^7(^9 ® %) with V g unitary. Now consider 
_ff 7+d7 with d7 I small, and expand its ground state 
as 

|M[7Vd 7 ]> = Vl - A 2 |m[?> 7 ]) + A|x 7 ) 
with (x 7 |/*[P 7 ]) = 0. With H 1+dl = H 1 + dH, 

= (A t [7 , 7 +d 7 ]|^ 7 +d 7 |/^[7 : ' 7 +d 7 ]) 
= A 2 (x 7 |ff 7 |x 7 ) + {n[V^}\dH\fi[V^}) (12) 
> A 2 A - ||dff|| , 

where A is the spectral gap of H 1 . Since A = ||/i["P 7 ]) — 
|/i.[7- > 7 +d 7 ]) I and dH — > as d7 — s- 0, this shows that 



1/ifP-y]) is continuous in 7, and since |/i[P 7 ]) is a poly- 
nomial in Vj, it follows that V 1 can be chosen to be a 
continuous function of 7 as well. 

Let us now study how the virtual representation of the 
symmetry is affected by a continuous change "P 7 +d 7 = 
V 1 + dV. Let us first consider the case where V 1 and 
dV are supported on the same virtual space, i.e., the 
bond dimension does not change, and let us restrict the 
discussion to the relevant space. The representation of 
the symmetry on the virtual level becomes 7 

z g ® z; = v-l Al u g v 1+Al , 

where Z g = Z g (X) is invertible, and Z* = (Zg 1 ) 1 ", 
cf. Eq. ([5]). Since both 'P 7 +d 7 and its inverse change con- 
tinuously, one can find a gauge such that this also holds 
for Z g . 

It remains to see that a continuous change of the repre- 
sentation Z g does not change its cohomology class; note 
that the choice of gauge for "P 7 +d 7 , Eq. ([5]), leads to a 
transformation Z g «-> Y Z g Y~ x which does not affect the 
cohomology class. Let us first assume that Z g = Z g (X) 
is diffcrcntiablc, and let 

Zg = Ug + U' g d 1 , (13) 

and U g Uh = e lui ~ 9 ^U gh . Then, we can start from 
Z g Z h = e iM9,h)+u'(g,h)d 7 ] Zgh 

(note that smoothness of Z g implies smoothness of 
exp[iu)(g,h)] = tr[Z g Z/ l ]/tr[Z g / l ] ) and substitute (|13p . 
Collecting all first-order terms in d7, we find that 

U g U' h + U' g U h = e l ^U' gh + iu,'(g, h)e l ^U gh . 
Left multiplication with U^Ug 1 = e'^^U^ yields 

and by taking the trace and using its cyclicity in the 
second term, we obtain 

Lj'{g,h) = -i{(}>g + 4> h -(i>g h ) (14) 

with (j> g = trlUg^-Ug]: This proves that differentiablc 
changes of Z g can never change the cohomology class 
of U). 

In case Z g = Z g (X) is continuous but not differentiablc, 
we can use a smoothing argument: For any e, we can find 
a differentiable Z g (X, e) such that ||^ 9 (A, e) — Z g (X)\\ < e, 
and define ui\ ie (g, h) via 

e i^A9,h) =tT [Z g (X,e)Z h (X,e)Z;i(X,e)} . 

This ui\^(g,h) (and in particular its real part) varies 
again according to (fT4"|) for any e and thus does not change 
its cohomology class, and since it is 0(e)-close to u!\(g, h), 
the same holds for the cohomology class of Z g (X). 
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In order to complete our proof, we also need to consider 
the case where V-y+d-y is supported on a larger space than 
V 1 . (The converse can be excluded by choosing d7 suffi- 
ciently small.) This can be done by considering the sym- 
metry on the smaller space, and is done in Appendix [C] 
Together, this continuity argument shows that we cannot 
change the cohomology class of the symmetry U g on the 
virtual level along a smooth gapped path H 7 and thus 
completes the classification of phases with unique ground 
states in the presence of symmetries. 



G. Phases under symmetries: Systems with 
symmetry breaking 

Having discussed systems with unique ground states, 
we will now turn our attention to systems with symmetry 
breaking, i.e., degenerate ground states, corresponding 
to non-injective MPS. Recall that in that case, the MPS 
projector V is supported on a "block-diagonal" space 

A 

n = C °° ® C ° a ■ ( 15 ) 

— ■ rLa 

Before starting, let us note that different form the 
injective case, symmetry broken systems can be invari- 
ant under non-trivial projective representations as well. 
However, we can always find a blocking k such that the 
symmetry of the symmetry on the blocked system is rep- 
resented linearly (see Sec. IIVC|) . and we will consider 
that scenario in the following. 



1. Induced representations and the structure of systems 
with symmetry breaking 

Let us first explain how the physical symmetry is real- 
ized on the virtual level — we will see that it has the form 
of a so-called induced representation of a projective repre- 
sentation (the proof can be found in Appendix |D|) . Con- 
sider an non-injective MPS l/x^]) and its parent Hamil- 
tonian. Then, any invariance of the Hamiltonian under 
a linear unitary representation of a group G can be un- 
derstood as invariance under an equivalent linear repre- 
sentation U g which — with the correct gauge for V, and 
in the correct basis — acts on the virtual system as 

a a£a 

Here, P g is a permutation representation of G permuting 
blocks with different a's in (fT5j) . P g leads to a natural 
partitioning of {1, . . . ,*4} into minimal subsets a invari- 
ant under the action of all P g which we call irreducible] 
the first direct sum in (|16p runs over those irreducible sets 
o. The V£ are unitaries, where h = h{g, a) is a function 
of g and a. Before explaining their algebraic structure, 



note that P g can be thought of as composed of permuta- 
tions Pg acting on irreducible subsets, and (fl6|) can be 
rewritten as a direct sum over irreducible subsets, 

^=0^(0^°®^) • (17) 

In the following, we will describe the structure of the 
symmetry for one irreducible subset a; in fact, degen- 
eracies corresponding to different subsets a are not of 
particular interest, as they are not protected by the sym- 
metry (which acts trivially between them) and are thus 
not stable under perturbations of the Hamiltonian. For 
a single irreducible subset a, the symmetry 

P g((&V£®V£) (18) 

has the following structure (which is known as an induced 
representation): Fix the permutation representation P° , 
pick an element a?o € a, and define a subgroup H C G as 

H = {g : 7r g (a ) = a } 

where n g is the action of P g on the sectors, ( a ) = 
P g H a . Further, fix a projective representation V* of the 
subgroup H . Then, the action of Vfi can be boosted to 
the full group G in (fT8|) by picking representatives hp 
of the disjoint cosets kpH, and labelling them such that 
Ttkp ((*o) = P- Then, for every g and a, there exist unique 
h and (3 such that 

gk a = kph , (19) 

and this is how h = h(g, a) in (|18|l is determined. Note 
that the action of the permutation is to map a to /3, so 
that (fT9|) carries the full information of how to boost the 
representation Vfi of the subset H to the full group; this 
is known as an induced representation. (It is straightfor- 
ward to check that this is well-defined.) 

Before classifying the phases, let us briefly comment on 
the structure of (fT7|) . The sectors a correspond to differ- 
ent symmetry broken ground states. The splitting into 
different irreducible blocks a corresponds to the break- 
ing of symmetries not contained in U g , and it is thus 
not stable under perturbations. Within each block, U g 
has subsymmetries which can be broken by the ground 
states — the a € a — and subsymmetries which are not 
broken by the ground states — corresponding to the sym- 
metry action Vfi <£> Vfi defined on subspaces H, where the 
symmetry acts as in the injective case. 

2. Structure of symmetry broken phases 

In the following, we will prove that two Hamiltonians 
with symmetry breaking are in the same phase iff i) the 
permutation representations P g are the same (up to re- 
labelling of blocks), and ii) for each irreducible subset a, 
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the projective representation Vfi has the same cohomol- 
ogy class. (In different words, we claim that two systems 
are in the same phase if U g permutes the symmetry bro- 
ken ground states in the same way, and if the effective 
action of the symmetry on each symmetry broken ground 
state satisfies the same condition as in the injective case.) 
Note that having the same P g allows us to even meaning- 
fully compare the Vfi as the subgroups H can be chosen 
equal. Also note that since the permutation is effectively 
encoded in the subgroup H, we can rephrase the above 
classification by saying that a phase is characterized by 
the choice of a subgroup H together with one of its co- 
homology classes. 

As a simple example, the Ising Hamiltonians H x = 
J2 a i a f+i an d H z = Y^ a i a i+i are m different phases 
if the Z2 symmetry (a x )® N is imposed: While H z can 
break the symmetry, H x can not and thus belongs to a 
different phase. 



will make use of a continuity argument. The argument 
will go in two parts: On the one hand, continuity implies 
that each of the symmetry broken ground states changes 
continuously, and thus the permutation action P g of U g 
stays the same. The effective action Vfi Vfi (modulo 
permutation) of the symmetry on each of the symmetry 
broken ground states, on the other hand, can be classified 
by reducing the problem to the injective case. 

Let us first show that the symmetry broken ground 
states change continuously. Let := \fi[P"]) be the 
(orthogonal) symmetry broken ground states of H 7 , with 
V" := V~{\-H a the restriction of P 7 to W a (this describes 
an injective MPS). For small changes d7, we can again 
(as in Sec. Ill F 3)) use continuity and gappedness of H 1 to 
show that the ground state subspaces of H y and H 1+ ^ l = 
H 1 + dH are close to each other, i.e., there exist ground 
states |x 7 +d 7 ) °f Hj+d-y such that 



3. Equality of phases 

The proof of equality requires again to devise an inter- 
polating path in the space of MPS which is sufficiently 
well-behaved (i.e., it yields a smooth and gapped path in 
the set of parent Hamiltonians) . It turns out that we can 
use essentially the same construction as in the injective 
case: We interpolate along a path 



M 7 )) 8JV © • • • © \u> A (l))® N = 1«. • • ■ »«>M7)> 



in a space with local dimensions J2 a (^a + £*q) 2 , where 



tt ( 7 )) = (i-7)ElM) + (i-7) £ M 

i=l i=D°+l 



the joint symmetry is given by 

p 9 ( 0(^ a © w h) ® © w£)) , 

a a£a 

where Vfi and W£ denote the projective representations 
for the two systems (as in the injective case, it can be 
thought of being embedded in the physical symmetry by 
blocking two sites). Again, the Hamiltonian along the 
path is the sum of the GHZ Hamiltonian and projectors 
of the form 1 — |w a (7))(u> a (7)| which couple to the GHZ 
degree of freedom; the resulting Hamiltonian is commut- 
ing and therefore gapped throughout the path. 



4- Separation of phases 

Let us now show that two phases which differ in ei- 
ther the permutation representation P g or the cohomol- 
ogy classes of the V£ cannot be transformed into each 
other along a gapped path. As in the injective case, we 



|IV> 7 )-|x 7 +d 7 )| <0*(dH) 



where 0*(dH) goes to zero as dH goes to zero. Since 
the |x 7 +d 7 ) can be expanded in terms of the \ip^ +d J}, 
they are MPS, \x% dl ) = \t*[Q^ +dl }}- Using continuity 
of the roots of polynomials, we infer that we can choose 
|2 7 +d 7 ~^7 I < 0*(dH). On the other hand, this implies 
that the Q 7+ d 7 are almost supported on H ai and thus 
|2 7 +d 7 - ^ 7 +d 7 l < 0*(dH). Together, this shows that 



r 7 



P 7 Q +d7 | <0*(dH) , 



i.e., the V", and thus the symmetry broken ground states 
|?/; 7 ), change continuously. Since P g describes the permu- 
tation action of the physical symmetry U g on the symme- 
try broken ground states, which is a discrete representa- 
tion, it follows that P g is independent of 7, i.e., it cannot 
be changed along a gapped path H~ ( . 

In a second step, we can now break the problem down 
to the injective scenario. To this end, do the following for 
each irreducible block of P g : Fix the «o used to define the 
subgroup H and its representation VJf, restrict the phys- 
ical symmetry U g to g G if, and consider the injective 
ground state |/x[7-' 7 ]) and the correspondingly restricted 
parent Hamiltonian (this can be done by adding local 
projectors restricting the system to the subspace given 
by do). Since U g , g £ H, leaves l/ip 3 ™]) invariant up to a 
phase, it is a symmetry of the restricted parent Hamilto- 
nian; it acts on the virtual level as Vfi Cg> Vft. The results 
for the injective case now imply that it is impossible to 
change the cohomology class of V£, thus completing the 
proof. 
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III. TWO DIMENSIONS 
A. Projected Entangled Pair States 

1. Definition 

Projected Entangled Pair States (PEPS) form the 
natural generalization of Matrix Product States to two 
dimensions^ For V : {C D )® i -> C d , the PEPS \n[V]) is 
obtained by placing maximally entangled pairs \u)d) on 
the links of a 2D lattice and applying V as in Fig. [5] As 
with MPS, PEPS can be redefined by blocking, which al- 
lows to obtain standard forms for V , discussed later on. 
Parent Hamiltonians for PEPS are constructed (as in ID) 
as sums of local terms which have the space supporting 
the 2x2 site reduced state as their kernel. 



2. Cases of interest 

As in ID, each PEPS has an isometric form to which it 
can be continuously deformed, yielding a continuous path 
of 7-deformed Hamiltonians along which the ground state 
degeneracy is preserved. There are three classes of PEPS 
which are of special interest. 

First, the injective case, where V is injective, 
and |/*[P]) is the unique ground state of its parent 
Hamiltonian^ 

Second, the block-diagonal case, where (ker V) 1 - = 

a =i^a, with H a = span{|i,j, k,l) : C«-i < i,j,k,l < 
£ a }; this corresponds again to GHZ- type states and 
Hamiltonians with ^4-fold degenerate ground states. 
These systems arc closely related to the ID non-injective 
case — they exhibit breaking of some local symmetry, and 
the ground state subspace is spanned by \fi [P\ 
Third, the case where the isometric form of V is 

v = ^Tv g ®v g ®w g ®w g (20) 

g 

(the ordering of the systems is top-down-left- right), with 
V g and W g unitary representations of a finite group G 
containing all irreps of G at least once; this scenario 
corresponds to systems where the ground state degen- 
eracy depends on the topology of the system, and which 




FIG. 6. PEPS are constructed analogously to MPS by apply- 
ing linear maps V to a 2D grid of maximally entangled states 
\ud). 



thus exhibit some form of topological orderj^ in partic- 
ular, for V g and W g the regular representation of G, the 
isometric form of these PEPS describes Kitacv's double 
model of the underlying groupJ^ All these three classes 
have parent Hamiltonians at the isometric point which 
are commuting and thus gapped. 

B. Gap in two dimensions 

1. Gap in the thermodynamic limit 

The major difference to the case of ID systems is that 
it is much more difficult to assess whether the parent 
Hamiltonian is gapped in the thermodynamic limit, and 
examples which become gapless at some finite deforma- 
tion < 7crit < 1 of the isometric form exist. For in- 
stance, the coherent state corresponding to the classical 
Ising models on a hexagonal lattice at the critical tem- 
perature has critical correlations and is thus gapless,— 
while it is injective and therefore its isometric form is 
gapped. In fact, this example illustrates that a smooth 
change in V can even lead to a non-local change in the 
PEPS \fi[P}). 

Fortunately, it turns out that in some environment of 
commuting Hamiltonians (and in particular in some envi- 
ronment of the three classes introduced above) , a spectral 
gap can be proven. To this end, let H = ^ hi, hi > 1 
with ground state energy A m i n (iJ) = 0, where the condi- 
tion 

hihj + hjh, > -|(1 - A) (hi + hj) (21) 

holds for some A > (here, each h% acts on 2 x 2 plaquc- 
ttes on a square lattice); in particular, this is the case for 
commuting Hamiltonians. Then, 

H 2 = h 2 + hihj + hh > AH , (22) 

(where the second and third sum run over overlapping 
and non-overlapping hi, hj, respectively), which implies 
that H has a spectral gap between and A, cf. Ref. [l|. 

As we show in detail in Appendix [EJ condition (f2l"j) 
is robust with respect to 7-dcformations of the Hamilto- 
nian. In particular, for any PEPS l/^T^]) with commut- 
ing parent Hamiltonian (such as the three cases presented 
above), it still holds for the parent of Q® N \^\P\) as long 
as A m in (Q)/A max (Q) ~ 0.967. Thus, while considering 
the isometric cases does not allow us to classify all Hamil- 
tonians as in ID, we can still do so for a non-trivial subset 
in the space of Hamiltonians. 

C. Classification of isometric PEPS without 
symmetries 

Let us now classify the three types of isometric PEPS 
introduced previously in the absence of symmetries; to- 
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gether with the results of the previous subsection, this 
will provide us with a classification of quantum phases in 
some environment of these cases. 

1. Systems with unique or GHZ-type ground state 

As in one dimension, any injective isometric V can be 
locally rotated to the scenario where V = 1, this is, it 
consists of maximally entangled pairs between adjacent 
sites. This entanglement can again be removed along a 
commuting path h e , Eq. ([5]), as in one dimension. This 
implies that any injective PEPS and its parent Hamilto- 
nian which is sufficiently close to being isometric is in the 
same phase as the product state. 

The case with block-diagonal V, such as for GHZ- like 
states, is also in complete analogy to one dimension: Up 
to a rotation, it is equivalent to the *4-fold degenerate 
GHZ state with additional maximally entangled pairs 
between adjacent sites, whose bond dimension D a can 
again couple to the (classical) value of the GHZ state. 
This local entanglement can again be removed along a 
commuting path, as in one dimension, and we find that 
all block-diagonal PEPS which are close enough to being 
isometric can be transformed into the ^4— fold degenerate 
GHZ state along a gapped adiabatic path. 

2. Systems with topological order 

What about the topological case of Eq. (|20|l ? Of 
course, additional local entanglement \ojd) can be present 
independently of the topological part of the state, which 
corresponds to replacing V g by V g ® 1 (and correspond- 
ingly for W g ), and this entanglement can be manipulated 
and removed along a commuting path. 

However, it turns out that the bond dimension D of the 
local entanglement can couple to the topological part of 
the state even though there is no local symmetry break- 
ing. In particular, the bond dimension D a can couple to 
the irreps R a (g) of V g and W g , i.e., we can change the 
multiplicity D a of individual irreps R a {g), 

0i? Q (.9) <->• 0^(5)®l Do • 

a a 

The interpolation between different multiplicities D a can 
be done within the set of commuting Hamiltonians, by 
observing that the Hamiltonian consists of two commut- 
ing parts)^ One ensures that the product of each irrep 
around a plaquette is the identity, and the other controls 
the relative weight of the different subspaces and thus 
allows to change multiplicities. The underlying idea can 
be understood most easily by considering a two-qubit toy 
model consisting of the two commuting terms 

h„ = \{t-Z®Z) , 
h x (6) = Af 2 1(1- X® X) Af 2 , 



where Ag = ( e a 9 ) : The term h z enforces the even-parity 
subspacc a|00) + /3|11), while h x (9) takes care that the 
relative weight within this subspace is |00)+6I 2 |11), which 
allows to smoothly interpolate between |00) and 1 00) + 
1 11) within the set of commuting Hamiltonians. 

Together, this proves that for a given group G, all 
PEPS of the form ([20|) . with representations V g and W g 
which contain all irreducible representations of G, yield 
PEPS which are in the same phase. On the other hand, 
it is not clear whether the converse holds: Given two fi- 
nite groups G, H with corresponding representations V g , 
W g and VI, W' hl for which Eq. (|2TJ)) yields the same (or 
a locally equivalent) map V — which means that the two 
models are in the same phase — is it true that the two 
groups are equal? While we cannot answer this question, 
let us remark that since both models can be connected 
by a gapped path, one can use quasi-adiabatic continu- 
ation^^ to show that their excitations need to have the 
same braiding statistics; this is, the representations of 
their doubles need to be isomorphic as braided tensor 
categories. Note that in Ref. the map V is used to 
map doubles to equivalent string-net models. 



3. More types of local entanglement 

Let us remark that while we have characterized the 
equivalence classes of isometric PEPS for the three afore- 
mentioned classes, this characterization is not complete, 
even beyond the difficulty of proving a gap: There are 
PEPS which can be transformed to those cases by local 
unitaries or low-depth local circuits, yet V has a differ- 
ent structure. The reason is that unlike in ID, local en- 
tanglement need not be bipartite. E.g., one could add 
four-partite GHZ states around plaquettes: while this 
is certainly locally equivalent to the original state, it will 
change the kernel of V, since only bipartite maximally en- 
tangled states can be described by a mapping V — > V®1. 
Thus, the previous classification can be extended to a 
much larger class of isometric tensors, by including all 
symmetries of ker V which can arise due to adding local 
entanglement. 



D. Symmetries in two dimensions 

How does the situation change when we impose sym- 
metries on the system and require the Hamiltonian path 
to commute with some unitary representation U g ? Sur- 
prisingly, imposing symmetries in two dimensions has a 
much weaker effect than in one dimension, as we will 
show in the following. In particular, we will demonstrate 
how to interpolate along a symmetry-preserving path be- 
tween arbitrary injective PEPS, and between any two 
GHZ-type PEPS given that the permutation action of 
the symmetry on the symmetry broken ground states is 
the same; note that the symmetry can in particular stabi- 
lize the degeneracy of GHZ-type states. Recall, however, 
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that in the following wc only show how to construct con- 
tinuous paths of PEPS; in order to turn this into a clas- 
sification of phases under symmetries, we need to restrict 
to the regions characterized in Sec. MI Bl where we can 
prove a gap. Yet, the following arguments show that the 
reasoning used for one-dimensional systems with symme- 
tries will not apply in two dimensions, and a more refined 
framework might be needed. 

1. Systems with unique ground states 

Let us start by studying the injective case. There, 
it has be showr*2i that any two maps V and V' which 
describe the same PEPS are related via a gauge trans- 
formation 

V 1 = V{e i<l 'Y ®Y* ®Z®Z*) , (23) 

with Y* = (F _1 ) T . This implies that any unitary invari- 
ance of l^tfP]) can be understood as a symmetry 

Ug = V~ X V g V = Vg®Vg®Wg®W g 

acting on the on the virtual system, with V g and W g 
projective unitary representations. 

While this is in complete analogy to the one- 
dimensional case, there is an essential difference: The 
representation of the symmetry on the virtual level is 
not invariant under blocking sites. By blocking k x £ 
sites, we obtain a new PEPS projector with symmetries 
V' = V® k and W' g = W® e , respectively. However, tak- 
ing tensor products changes the cohomology class, and 
in particular, for any finite-dimensional representation 
V g there exists a finite k such that V® k is in the trivial 
cohomology class. This is, by blocking a finite number of 
sites, any two PEPS |£t[Po]) an d \^[P{\) can be brought 
into a form where the symmetry on the virtual level is 
represented by V g ° , W® , and Vg , W g which are all in the 
trivial cohomology class. 

At this point, we can proceed as in the ID case 
and construct an interpolating path which preserves the 
symmetry using non-maximally entangled states |w(7)) 
[Eq. (|10p ] , now with joint symmetry 

(Vg O Vg) (g) (Vg (B^)® ^ ®W^<g> ^ ffi W^) . 

Note that the whole construction can be understood as 
the sequential application of two one-dimensional inter- 
polations (since at the isometric point the horizontal and 
vertical directions decouple) and thus, all arguments con- 
cerning technical points such as the embedding in the 
physical space can be directly transferred. 

2. Systems with local symmetry breaking 

The results of Ref. for the relation of two maps V 
and V', Eq. (|2"3"|) . can be readily generalized to relate 



different representations of a PEPS with local symmetry 
breaking, i.e., GHZ- type states. It follows that for any 
such PEPS, a symmetry can be understood at the virtual 
level as 

U a = P a (0 V£ ® n ® Wfi ® WZ) . 

a a£a 

where again the P g permutes the sectors, the o are min- 
imal subsets invariant under P g , and the Vfi and Wfi, 
together with P g , form an induced projective representa- 
tion. 

The permutation action P g of the symmetry is invari- 
ant under blocking. Thus, we can apply the same type of 
continuity argument as in ID to show that different per- 
mutations P g label different phases. (Since this argument 
proves that there cannot be a gapped path of Hamiltoni- 
ans on any finite chain connecting systems with different 
P g , it holds independently of whether we can prove a 
gap on the Hamiltonian along the path, as long as the 
initial and final systems are gapped.) On the other hand, 
the projective representations behave under blocking just 
as in the injective case: They map to tensor products of 
themselves, and thus, we can choose a blocking such that 
the V£ and W£ are all in the trivial cohomology class. If 
moreover the P g are equal, we can construct an interpo- 
lating path of PEPS just as in one dimension. 



IV. DISCUSSION 

In this section, we discuss various aspects which have 
been omitted in the previous sections. 



A. Matrix Product States 

1. Injectivity and translational invariance 

In order to obtain injectivity, it is necessary to block 
sites (say, k sites per block). Thus, the notion of locality 
changes: For instance, a two-local path of parent Hamil- 
tonians is 2fc-local on the unblocked system. Also, along 
such a path we can maintain translational invariance only 
under translations by k sites, i.e., on the blocked system. 

However, the number of sites which need to be blocked 
to obtain injectivity is fairly small, namely O(log-D) for 
typical cases^ and 0(D 4 ) in the worst case.— In partic- 
ular, this is much more favorable than what is obtained 
using renormalization methods^ where in addition to 
injectivity one has to go to block sizes beyond the corre- 
lation length. 

It should be noted that several of our arguments apply 
without blocking, thus strictly preserving translational 
invariance. In particular, the way in which symmetries 
act on the virtual level is the same without blocking, and 
the impossibility proofs for interpolating paths also ap- 
ply equally. On the other hand, it is not clear whether 
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we can construct interpolating paths without reaching 
injectivity of V: While we can move to isometric V along 
gapped paths, the structure of those isometric points is 
now much more rich; for instance, the AKLT projector 
is of that form. Also, certain alternative definitions of 
phases under symmetries (cf. Sec. IIVB[) behave differ- 
ently under translational invariance. Finally, in the case 
of symmetry broken systems, imposing translational in- 
variance can result in 'projective representations as phys- 
ical symmetries (this is discussed in Sec. HVCp . which 
leads to a much more involved structure of the ID in- 
duced representation in Eq. (|D9I) . 



2. Different MPS definition 
Typically, MPS are defined using a set of matrices A 1 

as 

]T tv[A i ^--A i -]\i u ...,i N ) . 

This definition can be easily related to our definition in 
terms of MPS projectors via 

Injectivity of V translates to the fact that the {A l }i 
span the whole space of matrices, and the "block- 
diagonal" support space in the non-injective case cor- 
responds to restricting the A 1 's to be block-diagonal and 
spanning the space of block-diagonal matrices. 

The effect of symmetries in this language can be writ- 
ten as follows: Eq. ([9]) becomes 

J2(U 9 )ijA j = V g T A% , 

i 

and Eq. ([16]) 

Y^ffig)^ = PlV?A%P g , 
3 

where the P g now permute the blocks of the A l, s, and 
together with V g form an induced representation. 

While the matrix formalism using the A l, s is more 
common, we choose the projector formulation since we 
believe it is more suitable for the purposes of this pa- 
per (with the exception of describing the block structure 
in the non-injective case): the V, or parts of it, are the 
maps which we use to conjugate the Hamiltonian with 
to get to the isometric form, and which we conjugate the 
U g with to obtain the effective action of the symmetry on 
the bonds. Also, in this formulation the isometric point is 
characterized simply by maximally entangled pairs with 
U ®U symmetry, rather than by an MPS with Kronccker 
delta tensors, and is thus more intuitive to deal with. 



B. Definition of phases 

1. Different definitions of phases 

There are other definitions of quantum phases: For 
instance, instead of gappedness one can ask for a path 
of Hamiltonians along which the ground states do not 
change abruptly. One can also right away consider 
the ground states instead of the Hamiltonian and ask 
whether two states can be transformed into each other us- 
ing (approximately) local transformations.— Both these 
definitions arc implied by ours: The existence of a gapped 
path implies that the ground states change smoothly, 
cf. Eq. (|12|). and using quasi-adiabatic continuation,— any 
gapped path of Hamiltonians yields a quasi-local trans- 
formation between the ground states^ 



2. Local dimension and ancillas 

In our definition of phases, we allowed to compare sys- 
tems with different local Hilbert space dimension. One 
way to think of this is to consider the smaller system as 
being embedded in the larger system. A more flexible 
way is to allow for the use of ancillas to extend the lo- 
cal Hilbert space. In fact, these ancillas are automatically 
obtained when blocking: Recall that we restricted our at- 
tention to the subspace actually used by V\ the remaining 
degrees of freedom (if sufficiently many to allow for a ten- 
sor product structure) can be used to construct ancillas. 
An explicit way to do so is to block two isometric tensors 
together: The state now contains a maximally entangled 
pair \ujd) (correlated to the GHZ in the non-injective 
case) which can be considered as a D 2 -dimcnsional an- 
cilla system. 

Note that we can, in the same way, obtain ancillas for 
systems with symmetries: After blocking three isometric 
sites, the maximally entangled states in the middle are in- 
variant under V g <S)V g <S)Vg<S)V g . Since also two maximally 
entanged states between sites (1,4) and (2,3) are invari- 
ant under that symmetry, the symmetry acts trivially on 
this two-dimensional subspace which thus constitutes an 
ancilla qubit not subject to the symmetry action. 



C. Symmetries 

1. Definition with restricted symmetry representations 

When defining phases under symmetries, we have 
allowed for arbitrary representations of the symmetry 
group along the path. What if we want to restrict to 
only the representation of the initial and final symme- 
try, Ug and Ug, respectively? It turns out that does 
not pose a restriction for compact groups, as long as 
at least one of the effective representations U g or U g 
after blocking to the normal form is faithful. Namely, 
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given such a faithful representation U g , we have that 

Xu(g) = tTUg = \tl[V g }\ 2 > (With Vg <g> Vg the Vff- 

tual realization of Ug), which implies that any represen- 
tation W g is contained as a subrepresentation in U® N 
for N large enough. (For finite groups, this follows as 
the multiplicity r^r J2 g Xu (g) N Xw(g) is dominated by 
Xw(l) since |xt/(ff)l is maximal for g = 1, cf. Ref. H3; 
for Lie groups, this argument needs to be combined with 
a continuity argument, cf. Ref. l28l) . Thus, starting from 
the symmetry representation Ug or Ug, we can effectively 
obtain any representation needed for the interpolating 
path by blocking N sites, proving equivalence of the two 
definitions. 



2. Definition with only one symmetry representation 

If the two systems to be compared are invariant un- 
der the same symmetry Ug = Ug, we might want to 
build a path invariant under that symmetry, instead of 
considering the symmetry U g = Ug © Ug as in our def- 
inition. In fact, these two definitions turn out the be 
equivalent, since we can easily map a path with symme- 
try U g to one with symmetry Ug = Ug, and vice versa. 
Given a path \ip(X)) (with corrsponding Hamiltonian) 
with symmetry Ug, we can add an unconstrained an- 
cilla qubit (cf. Sec. IIVB 2[) at each site and consider the 
path \ip(X)) ® (VT - A|0) +VA|l))® Ar , thus embedding 
the system in a space with symmetry U g . 

Conversely, we can map any path with symmetry U g to 
a path with symmetry U® = Ug. To this end, we can use 
Ug = Ug<E>1.2 to interpret the path as a path involving one 
system with [/^-symmetry and one unconstrained ancilla 
qubit per site, which again can be understood as being 
part of a (?7°) (g,2 -invariant subspace, cf. Sec. IIVB 2l Note 
that the factorization U g = U® ® I2 requires Ug and Ug 
to have the same phase, which motivates why we chose 
to lock the phase between Ug and Ug in such a scenario. 

3. Definitions with different classifications 

Defining phases in the presence of symmetries is suble, 
as different definitions can yield very different classifi- 
cations. In the following, we discuss some alternative 
definitions and their consequences. 

One possibility would be not to allow for an arbitrary 
gauge of the phases of the symmetries in Ug and Ug, 
but to keep the gauge fixed. In that case, the ID rep- 
resentations in (|B1|) and (|D1|) can enter the classifica- 
tion: While finite representations can be still removed 
by blocking sites, continuous representations remain dif- 
ferent even after blocking and cannot be changed into 
each other continuously, and thus, phases are addition- 
ally labelled by the continuous ID representations of the 
symmetry. 

This classification changes once more if one allows to 



compare blocks of different lengths for the two systems: 
E.g., for U(l) any two continuous representations e™ " 5 
and e mi ^ of $ € U(l) can be made equal by blocking \ni\ 
and | no | sites, respectively, as long as the signs of no and 
ni are equal. Therefore, in that case the phases are ad- 
ditionally labelled by the signs of the ID representations; 
this case has been considered in Ref. [Til 

To give an example where one might want to fix the 
phase relation between Ug and Ug , consider the two prod- 
uct states 1 0000 • • ■ ) and |0101 • • • ) under U(l) symme- 
try: Both Ug and Ug arise as subblocks of ® 
[with i?^ = exp(i(j>Z/2) the original U(l) representation], 
which suggests to fix the phase relation and in turn sepa- 
rates the two phases; note that this can be seen as a way 
to reinforce translational invariance. 



4- A definition with only one phase 

We are now going to present an alternative definition 
of phases under symmetries which yields a significantly 
different classification, namely that all systems with the 
same ground state degeneracy are in the same phase, just 
as without symmetries. The aim of this discussion is to 
point out that it is important to fix the representation of 
the symmetry group, and not only the symmetry group 
itself, in order to obtain a meaningful classification of 
phases in the presence of symmetries. 

We impose that along the path, the system is invariant 
under some (say, faithful) representation of the symmetry 
group, which may however change along the path. Then, 
however, it is possible to transform any state to pairs of 
maximally entangled states between adjacent sites in the 
same basis, and therefore to the same parent Hamilto- 
nian, by rotating the isometric form. Thus, all injective 
systems can be transformed to the same Hamiltonian pre- 
serving symmetries, and similarly for symmetry broken 
systems. It follows that even with symmetries, all sys- 
tems are in the same phase as long as they have the same 
ground state degeneracy. 

Note that while in our approach, we also use rotations 
to bring the system into a simple form, these "rotations" 
should just be thought of as choosing a convenient basis, 
and not as actual rotations. In particular, due to our 
way of imposing the two symmetry representations on 
orthogonal subspaces, U g = Ug © Ug, we never need to 
fix two different bases for the same subspace. 



5. Projective symmetries 

While our definition of phases under symmetries can 
be applied to both linear and projective representations, 
we found that in the injective case, symmetries are in 
fact always linear (cf. Appendix [Bj. In the non- injective 
case, however, there exist systems having projective sym- 
metries, such as the Majumdar- Ghosh model, which has 
the ground states (|01) — |10))' 8Ar and the same state 
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translated by one lattice site. This model is invariant 
under SU(2)/Z 2 [which can be understood as a projec- 
tive representation of SO(3)] and under the Pauli matri- 
ces (a projective representation of U>2 = Z2 x Z2). On 
the other hand, any d-dimensional projective represen- 
tation becomes linear (up to trivial phases) after taking 
its d'th power; thus, classifying phases under linear sym- 
metries as we did is in fact sufficient. Alternatively, any 
projective representation can be lifted to a linear repre- 
sentation, which is still a symmetry of the Hamiltonian: 
E.g., any J7 g -invariant Hamiltonian, [H,Uf N ] = 0, is 
also invariant under e*^U g , and thus under the represen- 
tation Vfc = k of the group K = (U g ) generated by the 
U g by itself. 

On the other hand, in case we want to build a joint 
symmetry representation U g = U g © U g before blocking, 
and do not want to lift the joint representation to a larger 
group, we get constraints on U g and U g : Firstly, both of 
them need to be in the same cohomology class — while 
the cohomology class can be changed by blocking, one 
might want to compare systems with a particular notion 
of locality. Also note that even if both symmetries are in 
the same cohomology class, one needs to adjust the trivial 
phases such that u(g,h) is actually equal up to a one- 
dimensional representation of the group, since otherwise 
U g does not form a representation; note however that 
this is actually a consequence of our requirement that 
U g forms a representation and does not follows from an 
underlying symmetry. 



6. Multiple copies and phases as a resource 

An interesting observation is that the classification of 
ID phases under symmetries is not stable if one takes 
multiple copies. For instance, one can construct a path 
of smooth gapped Hamiltonians which interpolates from 
two copies of the AKLT state to the trivial state while 
preserving SO (3) symmetry. More generally, for any two 
states |/4Po]) an ^ IM^i]) there exist k and fci such 
that |At[7'o])® feo can be converted to \fi[ Vi])® k l ■ This 
follows from the observation made in Sec. IIII Dl for two- 
dimensional systems: Taking tensor products changes the 
projective representation on the bond, and it is always 
possible to obtain a linear representation by taking a fi- 
nite number of copies. 

From a quantum information perspective, this shows 
that MPS which belong to different phases should not be 
regarded as a resource such as entanglement, but rather 
as characterized by conserved quantities such as parity 
(i.e., described by a finite group). The minimum require- 
ment for a resource should be that it cannot be created 
"for free" (e.g., by the quasi- local evolution created by 
a gapped path 9 -). However, an arbitrary even number of 
copies of the AKLT state can be created from one trivial 
state, which demonstrates that phases under symmetries 
should not be considered resources. 



7. Other symmetries 

While we have discussed the classification of phases for 
local symmetries Uf N , very similar ideas can be used to 
classify phases under global symmetries such as inver- 
sion or time reversal symmetry . 10 ' 11 ! 29 The fundamen- 
tal concept — that two V representing the same MPS or 
PEPS are related by a gauge transformation — equally 
applies in the case of global symmetries. However, it 
should be noted that there is an essential difference, in 
that the representation structure of the global symme- 
try need not lead to a representation structure on the 
virtual level, which in turn leads to classification crite- 
ria beyond cohomology classes. Let us illustrate this for 
reflection symmetry: Reflection is realized by applying 
a flip (swap) operator F to the virtual system, together 
with an operation 7r on the physical system reversing the 
ordering of the blocked sites. Thus, for an injective MPS 
l/xfT 3 ]) with reflection symmetry, we have that 

irVF = V{V-!®V-i) , 

where V-\ is the virtual representation of the non-trivial 
element of Z2 = {+1, —1}- (Note that if V is injective, 7r 
cannot be trivial, since otherwise F = V-i <g> V-i which is 
impossible; this shows that an injective MPS cannot have 
reflection symmetry unless it contains more than one site 
per block.) Applying a second reflection, we find that 

V = 7r(7r'PF)F = ttV (VI 1 <g> V-i)¥ 

= (irVF) (V- 1 ®V-i)=V(Y-iV-i®V-iV-i) 

i.e., the Z2 group structure of the symmetry is repre- 
sented on the virtual level as V-1V-1 = — similar 
to a projective representation, but with an additional 
complex conjugation. This relation allows for phases 
e i0 = j = i ; corresponding to symmetric and antisymmet- 
ric unitaries V—i, which cannot be connected continu- 
ously and thus label different phases; this observation 
has been used in Refs. [13 and to prove the separation 
of the AKLT phase from the trivial phase under either 
time reversal or inversion symmetry. 

D. Examples 

1. The six phases with D2 symmetry 

As an example for the classification of one-dimensional 
phases in the presence of symmetries, let us discuss the 
different phases under D2 = Z2 x Z2 symmetry, which 
appears e.g. as a subsymmetry of SO(3) invariant models; 
we will see that there is a total of six different phases 
under D2 symmetry. 

Let us label the elements of D2 = Z2 x Z2 by e = (0, 0), 
x = (1, 0), z = (0, 1), and y = (1, 1), with componentwise 
addition modulo 2. D2 = {e,x,?/, 2} forms a subgroup 
of SO(3) by identifying (1,0) with an x-rotation by tt 
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and (0, 1) with a z-rotation by w. There are two equiva- 
lence classes of projective representations, corresponding 
to the integer and half-integer representations of SO(3). 
In particular, the one-dimensional spin-0 representation 
pi = p x = Py = p\ = 1 belongs to the trivial class, 
and the two-dimensional spin-i representation p^, = X, 
Py = Y, and pi — Z (with X, Y, Z the Pauli matri- 
ces) belongs to the non-trivial class. For the following 
examples, we will always consider systems with physical 
spin S = 1; we label the basis elements by their S z spin 
component, | — 1), |0), and |1), and denote the physical 
representation of the symmetry group by 

R x = exp[inS x ] = - | - 1)<+1| - |0> <0| - | + , 
R z = exp[«rSJ = - | - 1)(-1| + |0)(0| - | + 1)(+1| , 
R y = cxp^nSy] = + | - 1)(+1| - |0)(0| + | + 1)(-1| . 

For systems with unique ground states, there are two 
possible phases. One contains the trivial state |0, . . . , 0), 
which can be trivially written as an MPS with bond di- 
mension D = 1, and 

Vt r w = |0) ; 

clearly, applying any physical transformation R w (w — 
x, y, z) to the physical spin translates to applying the one- 
dimensional representation 1 ® 1 on the virtual system. 
The corresponding Hamiltonian is 

ITtriv = -^|0}(0|+const . (24) 

i 

The second phase with unique ground state is illustrated 
by the AKLT stated which can be written as an MPS 
with D = 2, and a projector 

Vaklt =n s =i(l®*T) , 

with the projector onto the S = 1 subspace, and 

iY = (^J). Since we have that R x = H s =iX®XU s=1 , 
and correspondingly for y and z, it follows that the sym- 
metry operations are represented on the virtual level as 
R X V = T(X g> X), R y V = T(Y $ Y), and R Z V = 
V(Z<B)Z), which is a non-trivial projective representation 
of SO(3). The corresponding Hamiltonian is the AKLT 
Hamiltonian 

#AKLT = $ ' + 1^ ' ^+i) 2 + C ° nSt ] ■ 

i 

In order to classify all phases with symmetry breaking, 
we need to consider all proper subgroups of Z?2- There 
are four of them: 

H x = {e,x} 
H z ={e,z} 

Hy = { e >2/} and 
#triv = {e} 



The first three are isomorphic to Z2 , which has only triv- 
ial projective representations, and thus, each of them la- 
bels one phase; since H trlv also has only trivial represen- 
tations, it corresponds to a fourth phase with symme- 
try breaking. The number of symmetry broken ground 
states is \D2/H\, i.e., the first three cases have two-fold 
degenerate ground states, and the last case a four-fold 
degenerate one. 

Let us start with H z . A represcntant of that phase is 
a GHZ-type state of the form 

|GHZA = | + 1,...,+1) + |-L...,-1) , 

which can be written as an MPS with D = 2 and 

= | + 1)(0,0| + |-1)(1,1| , 

where the basis elements |0) and |1) correspond to two 
symmetry broken sectors. The action of the symmetry 
on the virtual level is 

R x V G hz z = -?ghz 3 [|0,0)(1,1| + |1,1)(0,0|] 
R,Vgbz. =-V G hz z [|0,0)(0,0| + |1,1)(1,1|] 

and correspondingly for R y ; i.e., while R z acts within the 
symmetry broken sectors, R x (and R y ) acts by permuting 
the different sectors. The corresponding Hamiltonian is 
the GHZ Hamiltonian 

i 

The same type of ground state and Hamiltonian is found 
for the subgroups H x and H y with correspondingly inter- 
changed roles (i.e., R x and R y , respectively, do not per- 
mute the ground states). Note that these three phases 
are indeed distinct in the presence of D2 symmetry, as 
there is no way how to smoothly change the element of 
the symmetry group which does not permute the sym- 
metry broken sectors. (This is related to the fact that 
D2I H is discrete, i.e., we are breaking a discrete symme- 
try; note that breaking of continuous symmetries does 
not fit the MPS framework since this would correspond 
to an infinite number of blocks in the MPS and would 
require gapless Hamiltonians.) 

Finally, choosing H tr i v = {e} gives a phase which fully 
breaks the D2 symmetry. A representative of this phase 
can be constructed by blocking two sites, with the four 
basis states 

\l) = \ + l)\+x) , |2) = | + l)|-x), 
|3) = |-l)|+x) ) |4} = |-l)|-*>, 

where | ± x) are the S x eigenstates with eigenvalues ±1, 

|±a?)oc|-l)±V2|0) + | + l) . 

We have that 

R x \±x) =-\±x) , R z \±l) = -|±1> , 
R x \±l) = -\Tl) , R z \±x) = -\ T x) . 
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It follows that all R w <S> R w (w = x, y, z) act as permuta- 
tions on the basis {|1), |2), |3), |4)}: 

{R X ®R X ) : |i) O |3); |2) o |4) ; 
(R Z ®R Z ) : |i) O |2); |3) O |4) ; 
(Ry®Ry): 1 1>**|4>; |2)o|3). 
Thus, the GHZ type state 

4 

|GHZ 4 > = X) I*.---.*) . 
fe=i 

which is an MPS with D = 4 and 

4 

7>GHZ 4 =53|fe)(fc,fc| , 
fc=l 

breaks all symmetries of D2, and represents yet another 
distinct phase with symmetry D2 = Z2 x Z2; the corre- 
sponding Hamiltonian is 

4 

#ghz 4 = -51 [5Z iMXMkj+i , 
j fe=i 

where the label j refers to the blocked sites. 

2. Phases under SO(3) and SU(2) symmetry 

Let us now consider symmetry under rotational invari- 
ance, imposed either as SO (3) or as SU(2) symmetry. We 
will find that under SO (3) symmetry, there are two pos- 
sible phases, represented by the spin-1 AKLT state and 
the trivial spin-0 state, respectively; on the other hand, 
we will show that if we impose SU(2) symmetry, there 
is only a single phase, as the AKLT state can be trans- 
formed into the trivial state keeping SU(2) symmetry. 

Let us start with SO (3) symmetry. In order to compare 
the AKLT state to the trivial spin-0 state, we need a rep- 
resentation of S0(3) which contains both the spin-1 and 
spin-0 representation; we will denote the representation 
as 

RniO) =exp[i0n-S]©l . 

While we could start with such a symmetry representa- 
tion right away, let us discuss how to obtain the same set- 
ting from a spin-1 chain by blocking: Blocking two spin-1 
sites gives a system with total spin 1 ® 1 = 2 © 1 © 0, con- 
taining both a spin-1 and spin-0 subspacc. It is straight- 
forward to check that after blocking two sites and apply- 
ing a rotation 1 ® iY, the isometric form of the AKLT 
state is 

7>aklt - 1 , (25) 

where the identity is on the two virtual spins with repre- 
sentation \® \ = 1 © 0. The physical rotation Ra acts 



on the virtual indices with the projective spin- ^ repre- 
sentation of the rotation group; it follows that the AKLT 
state is in the non-trivial equivalence class under SO (3) 
symmetry The trivial spin-1 state, on the other hand, 
is obtained by placing singlets between pairs of spin-1 
sites [i.e., between sites (1, 2), (3, 4), etc.]. After blocking 
these pairs, we obtain a product state with the spin-0 
state at each site, which is an MPS with D = 1 and the 
trivial projector 

Ariv = \S = 0) . (26) 

Thus, the rotation group Rfi{0) acts with the trivial rep- 
resentation on the virtual indices, and the trivial state is 
thus in a different phase than the AKLT chain. 

It should be noted that while D2 is a subgroup of 
SO(3), this does not imply that SO(3) exhibits all phases 
of D2 — indeed, the symmetry broken phases are miss- 
ing. The reason is that while for any subgroup H C D2, 
D 2 /H is finite, this is not true for SO (3). Since however 
SO(3)/H labels the symmetry broken ground states, this 
corresponds to breaking a continuous symmetry, which 
leads to gapless phases and cannot be described in the 
framework of Matrix Product States. 

Let us now turn our attention towards SU(2) sym- 
metry, and explicitly construct an interpolating path 
between the isometric projectors for the AKLT state, 
Eq. ([25]). and the trivial state, Eq. (|26]). Note that the 
difference is that now for both Paklt and Ptriv, the sym- 
metry action on the virtual level is a linear representation 
of SU(2) (namely the spin-i and the spin-0 representa- 
tion, respectively). We can now provide an interpolating 
path |* 7 ) = \lj{j))® n , with 

M 7 ))=7|0,0) + (l- 7 )(|l,l) + |2,2)) , 

where |\&o) corresponds to the isometric form ([25]) of the 
AKLT state, and to the isometric from (|26|) of the 
trivial state. Furthermore, for U <G SU(2), the whole 
path is invariant under the on-site symmetry 

(1®U)®(1®U) = 1®(U®U)®(U®U) (27) 

which contains the symmetries 1 and U ® U of the trivial 
and the AKLT state as subsymmetries; this proves that 
under SU(2) symmetry, the AKLT state and the trivial 
state are in the same phase. Note that the symmetry (|27j) 
is only a representation of SU(2), but not of SO(3), as 
integer and half-integer spin representations belong to in- 
equivalent classes of projective representations of SO(3); 
also, we cannot obtain this symmetry by starting only 
from the spin-0 and spin-1 representation of SU(2), as 
they are not faithful representations. Note that this in- 
terpolating path can already be ruled out by imposing 
a parity constraint on the total number of half-integer 
representations, by e.g. associating them to fcrmions. 
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V. CONCLUSION AND OUTLOOK 

In this paper, we have classified the possible phases 
of one-dimensional, and to a certain extent two- 
dimensional, systems in the framework of Matrix Product 
States and PEPS. We have done so by studying Hamil- 
tonians with exact MPS and PEPS ground states, and 
classifying under which conditions it is possible or impos- 
sible to connect two such Hamiltonians along a smooth 
and gapped path of local Hamiltonians. 

We have found that in the absence of symmetries, all 
systems are in the same phase, up to accidental ground 
state degeneracies. Imposing local symmetries leads to 
a more refined classification: For systems with unique 
ground states, different phases are labelled by equiva- 
lence classes of projective representations, this is, coho- 
mology classes of the group; for systems with degenerate 
ground states, we found that the symmetry action can 
be understood as composed of a permutation (permuting 
the symmetry broken ground states) and a representation 
of a subgroup (acting on the individual ground states), 
which together form an induced representation, and dif- 
ferent phases are labelled by the permutation action (this 
is, the subgroup) and the cohomology classes of the sub- 
group. In this classification, systems in the same phase 
can be connected along a path which is gapped even in 
the thermodynamic limit, while for systems in different 
phases, the gap along any interpolating path will close 
even for a finite chain. 

We have subsequently studied two-dimensional sys- 
tems and considered three classes of phases, namely prod- 
uct states, GHZ states, and topological models based on 
quantum doubles. We have shown that all of these phases 
are stable in some region, and demonstrated that within 
that region, and more generally within the framework 
used for MPS, imposing symmetries does not further con- 
strain the phase diagram. 

We have also compared different definitions of phases 
under symmetries and found that very different classifica- 
tions can be obtained depending on the definition chosen, 
ranging from scenarios where symmetries don't affect the 
classification at all, to scenarios where the classification 
is more fine-grained and e.g. the one-dimensional repre- 
sentations of the group partly or fully enter the classi- 
fication. In this context, it is interesting to note that 
there is a hierarchy in the classification of phases as the 
spatial dimension increases: Zero-dimensional phases are 
labelled by ID representations of the symmetry group 
(this is, its first cohomology group). This label vanishes 
in one dimension, and phases are now classified by the 
second cohomology group. This label, in turn, vanishes 
in three dimensions, and although we have demonstrated 
that we can't infer symmetry constraints from the conti- 
nuity of the PEPS projectors V alone, it is expected that 
phases under symmetries in two and more dimensions are 
still classified by higher order cohomology groups i 31 ' 32 

A central tool in our proofs has been the isometric form 
of an MPS or PEPS. Isometric MPS and PEPS are fixed 



points of renormalization transformations, and any MPS 
can be transformed into its isometric form along a gapped 
path in Hamiltonian space; this result allows us to restrict 
our classification of one-dimensional quantum phases to 
the case of isometric RG fixed points. Moreover, it gives 
us a tool to carry out renormalization transformations in 
a local fashion, this is, without actually having to block 
and renormalize the system; it thus provides a rigorous 
justification for the application of RG flows towards the 
classification of quantum phases. Let us add that the 
possibility to define an isometric form, as well as the 
possibility to interpolate towards it along a continuous 
path of parent Hamiltonians, still holds for not transla- 
tional invariant systems; however, without translational 
invariance we are lacking tools to assess the gappedncss 
of the Hamiltonian. 

Let us note that MPS have been previously applied 
to the classification of phases of one-dimensional quan- 
tum systems! 11 ' 29 ' 33 In particular, in Ref. |29j, MPS have 
been used to demonstrate the symmetry protection of the 
AKLT phase, and in Ref. [Til renormalization transforma- 
tions^ and their fixed points on MPS have been applied 
to the classification of quantum phases for one dimen- 
sional systems with unique ground states both with and 
without symmetries, giving a classification based on co- 
homology classes and ID representations. Beyond that, 
RG fixed points of PEPS have also been used towards the 
classification of phases for two-dimensional systems ! 14 ' 34 
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Appendix A: Gap proof for the ID path 

In the following, we show that the family of 7-deformed 
parent Hamiltonians which arise from the MPS path 
|/x['P 7 ]) interpolating between an MPS and its isomet- 
ric form is gapped. Recall that this family was defined 
as H 1 := Y^, h 7 (i, i + 1), with /i 7 := A 7 /ioA 7 > 0. 

We want to show that the path H 1 is uniformly 
gapped, i.e., there is a A > which lower bounds the 
gap of H-y uniformly in 7 and the systems size TV: This 
establishes that the |/i[P 7 ]), and the corresponding H~, 
are all in the same phase. To this end, we use a result 
of Nachtergaele^ (extending the result of Ref. [l] for the 
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injective case), where it is shown that any parent Hamil- 
tonian is gapped, and a lower bound on the gap (uniform 
in N) is given. 

In the following, we will use the results of Ref. 
to derive a uniform lower bound on the gap for all 
H y , 1 > 7 > 0. Let the MPS m atrices [Ai(i)] k i := 
Efc^W^lM)!^! (cf- Sec. HVA2ll : in the normal form, 
the Ai(j) have a block structure -A* (7) = © Af (7). Let 
E Q ( 7 ) := J2i A i(l) ® ^fh): and let |specE Q ( 7 )| = 
{A?( 7 ) > A£( 7 ) > ■■• > 0} be the ordered absolute 
value of the spectrum of £"(7) (not counting duplicates). 
Then, A 2 (7)^1(7) < 1, and since the spectrum is contin- 
uous in 7 6 [0; 1], and the degeneracy of Ai is A,— the ex- 
istence of a uniform upper bound 1 > r a > A 2 (7) /A" (7) 
follows. For a 7^ (3, let 



sup ■ 

X,Y 



(*[A a (>y);X]\*[Af> 
|$[A«(7);X])|||||$[^(7);y])| 



where \$[C;X\) := Eh,...,^ tr^ . . . C ip X]\h, . . . , i p ); 
i.e., ^^(7) is the maximal overlap of the p-site re- 
duced states of the MPS described by the blocks A" (7) 
and A»(i). With S a (i) := {E* *r[4f ( 7 )X]|i}|X}, and 
0(X,y) the maximal overlap between normalized vec- 
tors in the subspaces X and y, we have that Vt v a ^(7) < 

0(S Q (7)^,^(7)^) < 0(5 Q ( 7 ),5,3(7)) p . Moreover, 
since S a (0) _L 5^(0), and 5.(7) = Q 7 5.(0), we have 
that 



0(5„(7),5 /s ( 7 ))< sup 



MjgHl 

||<3 7 |v)|| ||Q 7 |iy)|| 



\M : 



12 



(u I ) =0 

V M11M22 



where M = irQ 2 ^ , it = 
submatrix of Q"i. For M > 0, 



(d + UXwL is some 2x2 



(M) ^ i A min (Q 7 ) 



< 1 



< 1 



M U M 22 " " A max (M) " A max (Q 7 ) 
< 1 - A mi „(Q 2 ) =: k < 1 , 

and we find that f2^ ^(7) < k p . Thus, there exists a p 
s.th. 



A{A - 1)k 
1-2(^-1) 



1 + D 2 tP 



and as Nachtergaele shows,- |A 2 p( 7 )(l — V^K P (^)) 2 is 
a lower bound on the spectral gap of H 1 . Here, A 2p (7) is 
the gap of H.y, restricted to 2p sites, which has a uniform 
lower bound as the restricted Hamiltonian is continuous 
in 7. This proves that has a uniform spectral gap for 
< 7 < 1. 



Appendix B: Standard form for injective MPS 
under symmetries 

In this section, we discuss how [/^-symmetry of an 
injective MPS is represented on the virtual level. To 
start with, it has been shown 7 that any two tensors V 
and V which (up to a phase) represent the same MPS 
can be related by a gauge transformation 

V' =V{e i ' t 'V ®V) . 

Given an MPS \n[P]) with t/ 9 -invariant parent Hamilto- 
nian (where U g is a linear or projective representation), 
it follows that |/^['P]) is invariant under U g up to a phase 
and thus, the action of U g on can be understood 

on the virtual level as 



or 



U g := V- 



%V = {e i ^V g ®V g ) 



(Bl) 



note that U g forms again a representation. From this, it 
follows that the V g form a projective representation, and 
in turn that <j) g forms a ID representation of G. This 
shows that U g , and thus U g , is a linear representation of 
G — systems with MPS ground states without symmetry 
breaking cannot be invariant under projective represen- 
tations. (Strictly speaking, we only find that <p g is in the 
same cohomology class as a linear representation, but in 
order to compare two systems, we need the same gauge 
for <j)g and <j) g , so we choose them to be linear represen- 
tations.) 

Under blocking, the physical symmetry U g in (|B1|) is 
mapped to Uf k , restricted to the range of the blocked 

map V® k \uiD)®^ k ~ 1 ' > ■ V g remains unchanged under block- 
ing; this suggests that it is suitable as a characteristic 
of a quantum phase (which we expect to be stable un- 
der blocking). The ID representation e 1 ^ 9 , on the other 
hand, changes under blocking to e lfc ^s . In particular, 
if the ID representation is finite, it can be removed by 
blocking, while this is not possible for continuous repre- 
sentations such as of U(l). 

However, our definition of quantum phases allows us 
to fully remove the ID representation e 4 * 9 in (|B1|) even if 
it is continuous, using the phase degree of freedom which 
we included in our definition. Consider first the scenario 
where we are free to choose the phase degree of freedom 
for the initial and final system independently. Then, we 
can choose to replace U g by e l ^ g U g , which will make the 
phase in (|B1|) vanish. Second, consider the case where 
the two physical symmetries are equal, U g = Ug =: U g . 
Then, we can choose the phase gauge such that the ID 
representation for one system vanishes, 

u° g = v g ° ® v g ° . 



On the other hand, since U„ 



U g , we have that 
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which also implies that e 1 ^ = 1 [e.g., by looking at any 
non-zero matrix element (i, i) x (J, j)]. Finally, if U g and 
Ug only share a subblock, this means that V® = X® © Y® 
s.th. Xg <£) Xg yields that subblock, and similary for Vg 1 , 
which again shows that e^n = 1. 



Appendix C: Continuity of cohomology class where 
subspace changes 

This appendix contains the proof omitted at the end 
of Sec. Ill F 31 (all the notation is the same as introduced 
there) : That the cohomology class obtained from "P 7 and 
V-y+d-f is the same even if "P 7 +d 7 is supported on a larger 
space than V T We will again show this for differentiable 
V-y, it extends to continuous by approximating them 
by a sequence of differentiable functions. Let Vy+d-y 

r . + v'd-f, with 



and 



V = 
I y 



P 




A B 
C D 



and note that the same block structure has to hold for 
the physical symmetry (as it is a symmetry for the whole 
path) , 



U„ 



Qg 
IL 



with R g , S g unitary. To first order in d 7 , the action of 
the symmetry on the virtual space is 



U 9 = V-l dl UgV 1+dl = T-'KgT + 0(d 7 2 ) 



with 
K„ 



T = 



(1 + X)Qg - XRg QgX - X Rg 

(1 + X)(R-Q) (t + X)Rg-QgX 



a o 

C- X D 



X = (P + Ad-f)- 1 BD~ 1 Cdj, and the "dressed represen- 
tations" 

Q g = (P + Ad 7 )~ 1 Q ff (P + Ad 7 ) , 
Rg = C 1 RgC . 

(This can be derived using the Schur complement to ex- 
press Vy+d-y an d can be readily checked by left multipli- 
cation with 7 3 7 _|_ ( j 7 .) It follows that the upper left block 
of K g K h is 

(1 + X)Q g Q h - XRgR h + 0(d 7 2 ) , 



i.e., the upper left block of U g (which corresponds to 
the virtual subspace used by 7 7 7 ) forms a representa- 
tion to first order in d 7 . The part corresponding to each 
of the two bonds thus forms a projective representation 
which changes smoothly in 7 , and which therefore cannot 
change its cohomology class. Since the cohomology class 
has to be the same for all subblocks, this completes the 
proof. 



Appendix D: Standard form for non-injective MPS 
under symmetries 

In the following, we will show that for any non-injective 
MPS l/xp 7 ]) with a parent Hamiltonian which is invari- 
ant under some local symmetry U g , the symmetry can 
be represented as the symmetry (fT6f of V . It is clear 
from ([2]) that any V with this symmetry will result in a 
Hamiltonian with the same symmetry. In the following, 
we will prove the converse. 

Using the same arguments as in Ref. l3ol one can show 
that any C/g-invariance of a parent Hamiltonian (at the 
moment, we are talking of a single unitary U g and only 
carry the subscript g for "future use" ) can be understood 
as resulting from an effective action 



u g :=r- 1 u g r = Pg 



0' 



P 9 V C 



(Dl) 



of the symmetry operation on the virtual system, using 
the correct gauge for V . Here, V" is unitary, the direct 



sum runs over the Hilbcrt spaces H a in (|15p . and P g per- 
mutes those Hilbcrt spaces. (Proof sketch, cf. Ref. l35l : 
C/ g -invariance of the Hamiltonian means that U g maps 
ground states to ground states. The different sectors 
Ron corresponding to different symmetry broken sectors, 
must be treated independently, since they don't interfere 
on any OBC interval. Thus, U g can on the one hand 
permute sectors and change their phase — this gives the 
P g and exp[z^g] — and it can on the other hand act non- 
trivially on each sector — since each of the sectors behaves 
like an injective MPS, this gives the V g a <g> V g a .) 

In the following, we will assume that P g acts irre- 
ducibly on the system in the sense that there are no sub- 
sets of {1, . . . ,„4} invariant under all P g . We can always 
achieve this situation by splitting (|D1|) into a direct sum 
over such irreducible cases. 

Let us now study what a linear representation struc- 
ture 



UgUh = Ug h 



(D2) 



(g, h G G) implies for the algebraic structure of P g , V g a , 
(We can always achieve a linear representation 



and 



by blocking; also note that the following argument can 
be generalized to projective representations.) First, since 
Pg is the only part of (|D1[) which is not block diagonal, 
it follows that the P g form a linear representation of G. 
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(Linearity follows since the entries of P g are and 1). 
Let us define 



sing that 7r fl (/3) = Kk y {nh(n k -i (J)))) = 7^(717,(0:0)) = 7- 



If we now substitute back 



vy g 



e*« V g a ® V g a 



(D3) 



Then, using (|D1[) the representation structure (|D2|) is 
equivalent to the relation 



w, 



(D4) 



for the Wg, where the permutations 717, arc defined 
via 'H- Kh [ a ) = Ph(Ha) and thus form a representation, 

Let us now show that (|D4j) implies that the U g can be 
understood as an induced representation; the following 
proof is due to S. Beigi^ Fix some ao, and let 



[Eq. dm])] in 



[Eq. (|D4|) ] . we find that forms a projective represen- 
tation of H, and e l< ^° a linear representation of H. 
Thus, we can write 



with 



tfU g K = U g D g 



H := {h : h £ G, Tr h (a a ) = a a } . 

Then, (|D4|) implies that Wf?° is a linear representation 
of H. We know we can write G as the disjoint union 
over cosets kpH labelled by the blocks ft = I,..., A, 
for a (non-unique) choice of kp £ G chosen such that 
Ttkpi&o) = ft- (This is where we need irreducibility.) We 
now have that 



" g 



^-i(/9) 



W: 



where we have used that n k -i(ft) = 7T fc = ao- Using 

the decomposition of G in cosets, we have that g and ft 
uniquely determine 7 and h £ H by virtue of 



gkp 

and thus, continuing (|D5I) . 



(D6) 



rh(Qo) vyr°w A/3 - 

7 ' 1 fe 



using that 717, (ao) = ao- We can now replace Wr_i using 
that 



1 



^ 1 ^wy iU>) w p 

kp 



Substituting this above, we finally obtain that 

. -1 



W 13 
vy g 



(D7) 



i.e., W^ 3 is fully determined by the representation W£° 
on H, together with the (arbitrary) unitaries W£° for all 
coset representatives kp. 
We now define a rotation 



then, in the rotated basis, U„ reads 



K%K = P g 



18 



18 



and 







(D8) 



(D9) 



where /;, = h(g,ft) is determined by (|D6j) . The diagonal 
operator D g acts independently on the different symme- 
try broken ground states of the system and thus com- 
mutes with the Hamiltonian; therefore, we can remove it 
by choosing the proper gauge for the physical symmetry 
according to our definition of phases under symmetries. 
If the symmetry of the initial and the final state overlap 
on a subsector of the ground state space, this implies (as 
for the injective case, Appendix [5} that the ID repre- 
sentations for this sector in D g are the same and thus 
can be removed by a joint gauge transformation. To- 
gether, this shows that for the classification of phases 
under symmetries in the non-injective case, any symme- 
try can be understood as a direct sum over independent 
sectors, where on each sector the action of the symme- 
try is given by (|D8j) . where h = h(g, ft) is determined by 
(ID61). 



Appendix E: Robustness of the 2D gap 

Here, we prove the robustness of a gap based on a 
condition of the form 

hihj + hjij > -f(l - + hf) ; (El) 

where we consider a square lattice with hi > 1 acting on 
2x2 plaquettes, Ay = A a for directly adjacent plaque- 
ttes i, j sharing two spins, and Ay = A^ for diagonally 
adjacent plaquettes i, j having one spin in common. (In 
Section IIII B| we have given the simplified version where 
A a = A d = A.) Then, 



H 2 = ~ h 1 + hihj + Y' hih ^ 



A a 



>hi 



<i3> 



>0 
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which implies a gap in the spectrum of H between and 
A = (A a + Ad)/2 > 0, and thus a lower bound A on the 
gap of H, cf. Ref. [J 

Let us now study the robustness of ()Elj) under 7- 
deformation of the Hamiltonian. Let hi, hj be projec- 
tors which satisfy hihj + hjhi > — i(l — Aij)(hi + ft?). 
(The proof can be modified for the /ii not being projec- 
tors.) Let hi be supported on systems AB, and hj on 
systems -BC, where the number of sites in systems A, B, 
and C is a, b, and c = a, respectively. (For the square 
lattice, a = b = c = 2 for directly neighboring terms, and 
a = c = 3, b = 1 for diagonally adjacent terms.) With 
Q 7 = (1 — 7)1 + 7Q < 1 as in the one-dimensional case, 



let A 7 ,x = (Q-y 1 ) with X = A,B,C and x = a,b,c. 
Then, the 7-deformed Hamiltonians are 

hi(j) = (A 7jj4 ® A 7iB )ft,i(A 7ij4 (g) A 7iB ) , 
^•(7) = (A 7) b (8 A 7) c)ft 3 -(A 7 ,s (8) A 7) c) ■ 

Let us define 

6 7 := A 7 B - 1 > , 
9 := A rain (Q) < 1 , 
^ 7 := ((l-7)+7A mi „(Q)r 2 > 1 , 

such that Q2q > an( j (,£ - 1)1 > 6 7 . Then, we 

find that 



hiA^shj+hjA^ghi + i/x" [l - &y + 8 (m 7 - 1)] (hi ® A- 2 C + A"^ ® hj) 

> hiA* >B hj + h 3 K 2 l B hi + § [1 - Ay + 8(/x 7 - 1)] (hi ® H c + Ha (8 

= /i i 6 7 /i J + ftj9 7 /ij + ftj [(/i 7 - 1)1] ftj + ftj [(/i 7 - 1)1] hj + h.hj + hjh, + |(1 - A y )(/i 1 + /i 3 ) 

> (h + h^e^hi + hj) > . 



By multiplying this with A 7jJ 4 ® A 7i b (g) A 7i c from both 
sides, we obtain a lower bound of type (|E1[) for the 7- 
deformcd Hamiltonian, 

^(7)^ (7) +hj(^hi (7) > -|[l-A ii ( 7 )][hi(7) + ^(7)] , 

(E2) 

with Ay (7) = ^^Ay- + (1 + 7/4 - 8^ +b ). This can be 
used to find an environment of any point in which the 
system is still gapped. In particular, in the case where 
the isometric parent Hamiltonian is commuting, and as- 



suming a square lattice, the lower bound on the spectral 
gap provided by (|E2[) is 

a, \ A a (7) + A d (7) 2 2 

A(7) = ^ = + M ^ + M7 - 2 AV • 

This gap vanishes at /i 7 w 1.07, limiting the maximal de- 
formation of the isometric tensor to A m i n (<5)/A max (Q) ~ 
0.967. 



1 M. Fannes, B. Nachtergaele, and R. F. Werner, Commun. 
Math. Phys. 144, 443 (1992). 

2 F. Verstraete and J. I. Cirac, (2004) , |cond-mat /0407066| 

3 F. Verstraete and J. I. Cirac, Phys. Rev. B 73, 094423 



(2006), cond-mat/0505140 
4 M 



Stat. Mech 



P08024 (2007), 
Rev. B 76, 035114 (2007), 



Hastings, 
larXiv:0705. 20241 

5 M. B. Hastings, Phys 
|cond-mat/0701055| 

6 B. Nachtergaele, Commun. Math. Phys. 175, 565 (1996). 

7 D. Perez-Garcia, F. Verstraete, M. M. Wolf, and J. I. Cir ac, 
Quant. Inf. Comput. 7, 401 (2007), |quant-ph/0608197| 

8 D. Perez-Garcia, F. Verstraete, J. I. Cirac, and M. M. W olf, 
Quantum Inf.Comp. 8, 0650 (2008). larXiv:0707.2260l 

9 M. B. Hastings and X. W en, Phys. Rev. B 72, 045141 
(2005) , |cond-mat/0503554| 

10 F. Pollmann, A. M. Turner, E. Berg, and M. Os hikawa, 
Phys. Rev. B 81, 064439 (2010). larXlv:0910.181ll 

11 X. Chen, Z. Gu, and X. Wen, Phys. Rev. B 83, 035107 



(2011). larXiv: 100 8.3745 

12 A. Kitaev, Ann. Phys. 303, 2 (2003), |quant-ph/9707u2l| 

13 F. Verstraete, D. Porras, and J. I. Cirac, Phys. Rev. Lett. 
93, 227205 (2004), [cond-mat/0404706l 

14 X. Che n, Z. Gu, and X Wen, Phys. Rev. B 82, 155138 
(2010) , larXiv: 1004.38351 

15 S. Michalakis and J. Pytel, (2011), [arXiv:1109.1588l 

16 S. Bravyi and M. B. Hastings, (2010). larXlv:1001.4363l 

17 F. Verstraete, J. Cirac, J. Latorre, E. Rico, and M. Wo lf, 
Phys. Rev. Lett. 94, 140601 (2005), |quant-ph/04102"27l 

18 C. H. Bennett, S. Popescu, D. Rohrlich, J. A. Smolin, 
and A. V. Thapliy al, Phys. Rev. A 63, 012307 (2000), 
|quant-ph /9908073 

19 N. Schuch, I. Cirac, and D. Perez-Garci'a, Annals of 
Physics 325, 2153 (2010). larXiv:1001.3807l 

20 F. Verstraete, M. M. Wolf, D. Perez- Garcia, and J. I. Cir ac, 
Phys. Rev. Lett. 96, 220601 (2006), |quant-ph/06010"75| 

21 M. B. Hasti ngs and T. Koma, Commun. Math. Phys. 265, 
781 (2006), math-ph/0507008 



24 



S. Bachmann, S. Michalakis, B. Nachtergaele, and R. Sims, 
(2011). larXiv: 1102.08421 

O. Buerschaper and M. Aguado, Phys. Rev. B 80, 155136 
(2009) , larXiv:0907.2670l 

D. Perez-Garcia, M. Sanz, C. E. Gonzalez-Guillen, M. M. 
Wolf, and J. I. Cirac, New J. Phys. 12, 025010 (2010), 
arXiv.org:0908.1674. 

D. Perez-Garci'a, M. Sanz, A. Cadarso, M. Wolf, and 
J. Cirac, in preparation . 

M. Sanz, D. Perez-Garcia, M. M. Wolf, and J. I. Cirac, 
IEEE Trans. Inf. Theory 56, 4668 (20091 . [arXTv:0909. 53471 
W. Fulton and J. Harris, Representation theory: a first 
course (Springer, 1991). 



http: / /mat hoverflow.net /questions/58633 

F. Pollmann, E. Berg, A. M. Turner, and M. Oshikawa, 

(2009) , arXiv.org:0909.4059. 

A. Affleck, T. Kennedy, E. H. Lieb, and H. Tasaki, Com- 
mun. Math. Phys. 115, 477 (1988). 

R. Dijkgraaf and E. Witten, Commun. Math. Phys. 129, 
393 (1990). 

A. Kitaev, private communication. 

L. Fid kowski and A. K itaev, Phys. Rev. B 83, 075103 

(2010) . larXTv:1008.4138l 

Z. Gu, Z. Wang, and X. Wen, (2010). larXlv:1010.1517l 
M. Sanz, M. M. Wolf, D. Perez-Garcia, and J. I. Cirac, 
Phys. Rev. A 79, 042308 (2009). larXlv:0901. 22231 

S. Beigi, private communication. 



